UNCLASSIFIED 


_ AD  NUMBER _ 

AD224373 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited.  Document  partially  illegible. 


FROM: 

Distribution  authorized  to  U.S.  Gov’t,  agencies 
and  their  contractors; 

Administrative /Operational  Use;  MAR  1955.  Other 
requests  shall  be  referred  to  Defense  Public 
Affairs  Office,  DFOISR,  Washington,  DC  20301. 
Document  partially  illegible . 


_ AUTHORITY _ 

rand  corp  ltr,  via  DoD  dtd  23  aug  1967 


THIS  PAGE  IS  UNCLASSIFIED 


UNCLASSIFIED 

ad  2  2  4  3  7  3 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION,  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  fomilated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


AN  OPTIMAL  INVENTORY  POLICY  FOR  A 
MILITARY  ORGANIZATION 

Edward  B.  Berman 
Andrew  J.  Clark. 

P-647 

March  30>  1955 


-7 6c  hx“H  I  I  l  J  (?<yijkoiati6H 


, P-647 

ft 

AN  OPTIMAL  INVENTORY  POLICY  FOR  A  MILITARY  ORGANIZATION 

Edward  B.  Berman  . 

Andrew  J.  Clark  =J 

"  1.  INTRODUCTION  • 


The  purpose  of  this  paper  is  to  present  what  is  believed  to  be  a 
realistic  optimal  inventory  policy  for  a  military  organization.  Our  prime 
consideration  has  been  practicability  of  application  rather  than  elegance  or 
generality.  To  further  this  end,  we  use  the  simplest  mathematical  techniques 
which  provide  the  desired  results  and  avoid  unnecessary  mathematical  rigour. 
If  we  overelaborate  in  describing  the  concepts  involved,  It  is  only  because 
we  wish  to  engender  more  complete  understanding  and  possible  application  of 
the  theory  in  the  "real  world". 

The  supply  system  which  we  consider  has  two .echelons  of  supply,  which  we 
shall  call  bases  and  depots,  and  are  not  unlike  the  wholesale  and  retail  out¬ 
lets  of  a  large  integrated  company.  Our  prime  assumption  for  the  base  stock- 
age  policy  is  that  the  concept  of  reorder  level  and  fixed  reorder  amount  is 
used.  Thus,  when  the  stock  oh  hand  at  a  base  for  a  given  item  reaches  the 

2/ 

reorder  level,  R,  an  order  is  submitted  to  a  depot  for  a  fixed  amount,  A 
Our  first  problem  will  be  to  find  the  R  and  A  which  allows  the  base,  as 
an  isolated  supply  facility,  to  operate  at  least  cost.  We  include,  both  the 
case  of  recoverable  items  (items  which  can  be  reclaimed  and  repaired)  and 


1/  While  the  authors  are  employed  by  The  RAND  Corporation,  Santa  Monica,  the 
views  expressed  in  this  paper  are  their  own. 

?/  For  conditions  under  which  this  policy  is  optimal,  see  On  the  Optimal 
Character  of  the  (s.S)  Policy  in  Inventory  Theory,  by  A.  Dvoretzky,  J.  Kiefer, 
and  J.  Wolfbwitz  in  Econometrica,  Vol.  ?1,  No.  1,  October,  1953,  pp  586-596. 
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non-recoverable  items.  As  a  matter  of  general  interest,  we  extend  our  re¬ 
sults  to  include  the  case  of  a  commercial  firm  using  the  reorder  level  con¬ 
cept.  In  our  base  results,  we  include  all  of  the  usual  costs  incurred  by 
the  base,  and  we  also  include  two  distinct  kinds  of  pipeline  time  (the  time 
between  the  issuance  of  a  requisition  and  the  receipt  of  the  materiel).  Our 
results  also  include  the  case  of  unknown  demand,  where  the  parameter  of  the 
demand  density  function  is  known  only  to  within  a  probability  function. 

We  then  extend  our  results  to  include  a  system  of  military  bases  and 
depots.  Here  we  discuss  three  general  classes  of  procurement  policies  which 
are  now  in  common  use  in  military  supply  systems,  These  are  life-of-type 
procurement,  wherein  sufficient  spare  parts  for  an  end  item  arc  purchased  at 
once  to  last  the  life  of  the  end  item;  periodic  procurement,  where  purchases 
are  made  at  regular  time  intervals  in  sufficient  quantity  to  satisfy  antici¬ 
pated  needs  during  the  period;  and  open-end  contracting,  where  an  order  on 
a  manufacturer  is  placed  each  time  the  system-wide  stocks  reach  a  certain 
minimum  level.  These  procurement  policies  are  characteristic  of  military 
organizations  due  to  the  allotment  of  funds  on  a  fisca.l  year  basis.  For 
each  procurement  policy,  we  include  the  case  of  recoverable  and  non-recover- 
able  items. 

Thera  are  two  additional  system  problems  for  which  no  solution  is. pre¬ 
sented.  These  are  the  problem  of  finding  the  optimal  ordering  period  for 
periodic  procurement,  and  the  problem  of  determining  optimal  delivery  periods 
after  procurement.  This  last  problem  is  represented  in  the  case  where  pro¬ 
curement  is  initiated  at  the  first  of  the  year  for  the  entire  year,  and  de¬ 
liveries  are  made  at  specified  times  and  in  specified  amounts  throughout  the 
year.  This  kind  of  procurement  and  delivery  policy  is  in  very  common  use. 

As  a  method  of  implementing  the  proposed  policies,  one  might  think  of 
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the  existence  of  an  information  processing  center  in  the  supply  system.  This 
information  processing  center  would  maintain  current  balances  of  each  item  on 
hand  at  each  of  the  bases  and  depots.  As  issues  are  made  on  the  bases,  the 
information  processing  center  is  so  notified,  and  balances  are  adjusted  at 
the  end  of  the  day  for  all  issues  during  the  day.  The  information  processing 
center,  having  calculated  the  reorder  level  and  amount  of  order  for  each  item 
on  each  base,  checks  each  balance  after  posting  issues  and  determines  whether 
or  not  the  balance  has  gone  below  the  reorder  level.  If  so,  the  information 
processing  center  issues  an  order  on  a  shipping  activity  to  ship  an  amount 
to  the  base.  The  time  between  the  issue  at  the  base  which  caused  the  balance 
to  fall  below  the  reorder  level  and  the  acceptance  of  the  materiel  at  the 
base  is  recognized  as  the  routine  pipeline  time.  If  a  base  should  exhibit  a 
need  for  an  amount  not  available  at  the  base,  (due  perhaps  to  an  unantici¬ 
pated  high  demand)  such  need  is  relayed  by  the  information  processing  center 
to  a  shipping  activity.  Tne  time  between  the  first  recognition  of  the  need 
at  the  base  and  the  actual  receipt  of  the  materiel  is  the  priority  pipeline 
time.  The  information  processing  center,  knowing  the  balances  available 
and  the  consumption  history  of  all  the  components  of  the  system,  is  in  a 
position  to  implement  the  policies  which  are  subsequently  derived.  The 
policies  are  not  contingent  upon  the  existence  of  such  an  information 
processing  center;  however,  the  operation  of  the  center  would  be  consistent 
with  the  policies. 

All  of  our  results  depend  upon  the  existence  of  demand  probability 
functions.  These  functions  may  be  difficult  to  obtain  in  practice.  They  may 
depend  upon  past  consumption,  future  programs,  operational  characteristics, 
and  numerous  other  factors.  On  the  other  hand,  any  requirements  and  levels 
computational  method  requires  such  functions  in  one  form  or  another.  Such 
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functions  are  now  being  used  whether  realized  or  not.  For  example,  the  simple 
procedure  of  estimating  the  next  sixty  day*s  demand  to  equal  the  past  sixty 
dayfs  consumption  automatically  establishes  a  probability  demand  function, 
even  though  rather  trivial  in  this  case.  Accepting,  then,  the  fact  that 
such  functions  must  be  established  and  used  for  any  method,  we  address  our¬ 
selves  to  the  problem  of  using  such  functions  to  provide  supply  support  in 
some  ’'optimal"  fashion.  We  assume  these  demand  probability  functions  to  be 
given;  their  actual  derivation  is  outside  the  scope  of  this  paper. 

The  general  technique  used  in  this  paper  establishes  decisions  based  on 
cost  considerations  alone.  We  assume  that  any  non-pecuniary  elements  in¬ 
volved  in  the  decisions  -  principally  the  end  item  being  out  of  commission 
for  parts  -  can  be  conveited  into  a  cost  factor.  Behind  this  assumption  lies 
a  further  assumption  to  the  effect  that  the  military  organization  can  purchase- 
more  end  items  as  an  alternative  to  more  logistics  support,  and  vice  versa. 
Thus,  given  a  programmed  requirement  for  the  number  of.  in-commission  end  items 
and  for  the  activity  of  these  end  items,  the  assumed  military  objective  is  to 
minimize  the  sum  of  the  costs  of  logistics  support  and  the  costs  of  a  pool 
of  end  items  being  out  of  commission  for  parts.  Alternatively,  given  a  pro¬ 
gram  of  end  item  activity  and  a  single  fund  limit  for  both  the  purchase  of 
whole  end  items  and  for  logistics  support,  the  assumed  military  objective  may 
be  expressed  as  the  maximization  of  the  number  of  end  items  in  commission. 

Using  either  statement  of  the  military  objective,  a  balance  may  be  preserved 
between  whole  end  items  and  logistics  support  by  charging  a  depletion  penalty 
against  the  logistics  decisions.  This  depletion  penalty  consists  of  the  total 
cost  of  the  end  item,  including  purchase  and  carrying  costs,  divided  by  the 
number  of  days, in  its  expected  useful  life,  and  is  assessed  for  each  day  the 
end  item  is  out  of  commission  for  parts.  This  depletion  Charge,  assessed 
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for  each  spare  part  which  could  cause  the  end  item  to  be  out  of  commission  „ 
for  lack  of  the  part,  is,  in  a  sense,  the  contribution  of  the  spare  part  to 
the  pool  of  end  items  out  of  commission*  Thug,  no  matter  how  many  of  the 
spare  parts  are  stocked,  there  is  some  probability  of  these  not  being  enough 
to  keep  the  end  item  in  commission,  and  therefore,  each  spare  part  has  some 
expected  contribution  to  the, pool  of  end  items  out  of  commission  for  parts. 
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PART  1  -  THE  BASE 

2 .  THE, GEHERAL  SOLUTION  FOR  THE  BASE1 

In  this  section  we  consider  the  base  as  separate  from  the  system*  Im¬ 
plicit  in  the  derivation  of  the  base  policy,  are  the  following  assumptions, 
most  of  which  have  been  already  suggested  in  the  introduction: 

A.  The  total  issues  which  occur  during  the  day  are  assumed  to  have 
occurred  one  at  a  time,  and  evenly  spaced  throughout  the  day. 

B.  The  balance  on  hand  is  compared  with  the  reorder  level  after  each 
issue,  and  a  requisition  initiated  if  the  balance  on  hand  is  equal  to  the  re¬ 
order  level.  The  requisition,  each  time,  is  for  an  amount,  a  ,  less  any 
amounts  received  on  priority  requisitions  during  the  previous  process  of  re¬ 
plenishment.  Also,  the  stock  on  hand  when  the  order  arrives  is  not  less  than 
the  reorder  level. 

Assumption  A  and  B  taken  together  permit  the  further  assumption  that 
the  balance  on  hand  is  exactly  equal  to  the  reorder  level  at  the  time  the 
requisition  is  submitted. 

C.  The  ability  of  depots  to  fulfill  requisitions  is  assumed  infinite. 
This  does  not  imply  that  such  is  the  case,  but  rather  that  any  costs  arising 
from  depot  failure  should  be  assessed  against  the  inventory  policy  of  the 
system  as  a  whole. 

Assumption  C  allows  us  to  look  at  the  base  as  an  isolated  supply 


The  general  approach  taken  in  the  base  solution  is  similar  to  that  of 
T.  M.  Whitin,  Theory  of. Inventory  Management,  pp  56-62.  We  are  indebted  to 
T.  M.  Whitin  for  a  critical  review  of  a  preliminary  draft  of  this  part  of 
the  paper.  We  are  also  indebted  to  R.  Bellman,  0.  Morgenstern,  and  others 
for  their  comments  on  the  preliminary  draft  of  this  paper.  Perhaps  we 
should  also  express  our  appreciation  to  K.  J.  Arrow,  T.  Harris,  and  J. 

Marschak,  for  it  was  their  paper,  optimal  Inventory  Policy.  Econometrica, 

Vol.  19,  1951,  duly,  pp  250-27?,  which  initially  introduced  this  problem  to  us. 
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activity  for  our  initial  derivation*  This  assumption  is  modified  when  we 
consider  the  base  as  a  component  of  the  system. 

D*  Whenever  demand  exceeds  the  balance  on  hand,  an  emergency  requisition 
is  submitted  for  the  excess  of  amounts  demanded  over  amounts  on  hand.  The 
emergency  requisition  requires  premium  communication  and  transportation. 

E.  The  routine  pipeline  time,  p,  and  the  emergency  pipeline  time,  p5 
are  assumed  to  be  constants. 

Next,  we  will  define  the  various  costs  encountered  in  the  operation  of 
the  base,  together  with  the  symbols  to  be  used  in  the  derivation  of  the  base 
policy.  First,  however,  we  should  clarify  our  position  in  one  respect.  We 
derive  our  stockage  policy  on  an  individual  item  of  supply  basis;  hence,  our 
various  costs  are  also  on  that  basis.  Since  many  of  the  costs  are  reduced, 
per  unit,  due  to  aggregation  of  items  (such  as  the  fixed  paperwork  cost  of 
requisitioning,  wherein  if  several  items  are  included  in  the  requisition,  the 
per  unit  cost  is  less  than  the  cost  of  a  single  item  on  the  requisition),  we 
assume  that  such  reduced  per  unit  costs  are  used,  In  practice ,  these  reduced 
costs  due  to  aggregation  can  be  found  by  sampling  techniques. 

The  symbols  and  component  costs  of  supply  used  in  subsequent  calcula¬ 
tions  are: 

(l)  k  «  fixed  costs  of  handling  the  paperwork  and  communications  for  a 
routine  requisition. 

(?)  A  *  fixed  amount  of  reorder. 

(3)  p^  +  PpA  a  packaging,  inspection,  and  handling  costs,  where  p^  is  the 

fixed  component  and  p  is  the  per  unit  cost  associated  with 
the  sise  of  the  order,  A. 

(4)  tj  *  tpA  a  cost  of  transporting  an  amount  A  from  the  supply  depot  to 


the  base. 
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(5)  y  58  amount  of  stock -on  hand. 

(6)  d-^  +  dpy  -  warehousing,  depreciation,  and  obsolescence  costs  per  unit 

of  time  for  holding  an  amount  y,  These  costs  are  also 
referred  to  as  "holding”  costs. 

(7)  o  =  amount  of  depletion. 

(8)  p  ®  routine  pipeline  time  (time  from  when  the  balance  on  hand  reaches 

the  reorder  level  to  the  time  of  receipt  of  the  materiel)  “  constant, 

(9)  p  58  priority  pipeline  time  (time  from  the  indication  of  need  for  an  item 

unavailable  at  the  base  to  the  time  of  receipt  of  the  item)  “  constant. 

(10)  bQ  +  b]_q  +  bppq  »  depletion  penalty,  or  the  cost  of  understocking  the  item. 

In  this  cost,  bQ  is  the  fixed  cost,  if  any,  which  is 
independent  of  the  amount  of  depletion.  The  cost  b-^ 
is  the  cost  associated  with  the  amount  of  depletion 
but  not  length  of  depletion.  The  component  bp p  q  is 
the  cost  incurred  through  loss  of  utility  of  an  end 
item  being  out  of  commission  for  p  days.  Here  b^ 
could  be  the  cost  per  day  of  having  an  extra  end  item 
available  for  use.  The  component  costs  of  the  deple¬ 
tion  penalty  will  be  discussed  in  further  detail  later. 

(11 )  x  «  number  of  the  item  demanded  per  unit  of  time.  In  subsequent  work, 

x  is  a  random  variable. 

(l?)  f(x)  *  demand  probability  (density)  function  of  the  random  variable,  x. 

This  function  may  be  either  continuous  or  be  defined  for  only 
integral  values  of  x. 

rx 

(13)  F(x)  B  I  f (t )dt  »  the  integral  or  cumulative  demand  probability 
''o 

function.  Since  f(x)  may  be  either  a  continuous' 
function  or  a  density  function,  such  as  a  Poisson 
distribution,  f(x)  must  be  integrable  in  the 


Stieltjes  sense  over  all  intervals  on  the 
positive  x-axis  with  zero  as  a  lower  limit*  We 
also  have  the  restriction  that  F(x)  1  as 

x  — ^  oo. 

(14)  R  =  reorder  level.  This  is  the  point  in  the  stock  inventory  where  a 

routine  requisition  is  required. 

(15)  0  3  time  period  between  the  arrival  of  two  successive  orders. 

We  notice  that  most  of  the  costs  encountered  by  the  base  (in  particular, 
items  (3),  (4)#  (6),  and  (10) above)  are  expressed  as  linear,  non-homogeneous 
functions  of  the  form  y  »  a  +  bx.  These  costs,  in  reality,  do  not  assume 
this  form,  but  may  be  reasonably  approximated  by  such  functions.  Actually, 
in  our  later  formulation,  we  could  just  as  easily  consider  these  costs  to  be 
expressed  as  arbitrary  functions  expanded  in  power  series.  However,  in 
practice,  probably  the  best  we  can  do  is  to  find  the  linear  approximations 
to  the  cost  functions;  hence,  we  restrict  ourselves  to  such  linear  cost 
functions  in  our  formulation. 

In  order  to  obtain  a  clearer  understanding  of  the  situation  and  to  . 
facilitate  the  calculation  problem,  we  shall  take  the  routine  pipeline  time 
as  the  fundamental  time  unit.  This  practice,  in  fact,  is  one  of  the  main 
characteristics  of  our  approach  and  provides  a  significant  simplification  of 
the  problem.  This  unit,  then,  will  also  serve  in  the  definition  of  those 
parameters  defined  in  terms  of  units  of  time,  such  as  the  parameters  d^  and 
d?f  the  variable  x,  and  the  requisitioning  period,  0.  Thus,  for  example, 
we  speak  of  9  as  being  "so  many  routine  pipeline  times". 

Our  base  inventory  as  a  function  of  time  may  now  be  portrayed  graphically' 
as  shown  in  Figure  (l).  Figure  (l)  represents  a  case  of  reality,  where  issues 
can  be  male  several  at  a  time  and  at  any  time  during  the  day.  Ry  applying 
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assumption  A,  we  amend  Figure  (l)  to  appear  as  shown  in  Figure  (2).  Thus, 
the  issue  of  three  items,  causing  a  drop  in  the  inventory  from  A  to  B  in 
Figure  (1),  is  represented  in  Figure  (2)  a3  three  separate  issues  of  one 
item  each  and  at  equal  intervals  throughout  the  day.  Also,  note  that  in 
Figure  (2),  the  routine  pipeline  time,  p  ,  is  considered  to  he  an  integral 
number  of  days;  this  is  not  at  all  necessary,  but  in  practice  such  would 
probably  be  the  case. 

We  are  now  in  a  position  to  write  a  function  representing  the  cost  of 
operation  at  the  base  for  a  given  item  and  for  a  typical  requisitioning 
period,  9; 


(1)  L(9,R,A)  -  (dp  +  dpy)  9  +  k  +  (pp  +  p?A)  +  (tp  +  t2A)  +  bQ  |jL  -  F(R)~| 

roo 

+  (b-^  +  b?p)  J  (x  -  R)  dF(x). 

R  , 

The  first  term  of  the  right  member  represents  the  cost  of  holding  a 

quantity  y  in  stock  over  a  period  of  time  9.  The  next  three  terms  are  the 

handling,  packaging,  and  transportation  costs  for  the  routine  requisition.  The 

term  bQ  |jL  -  F(R)^|  represents  the  expected  value  of  those  costs  of  depletion 

which  do  not  depend  on  the  length  of  depletion  nor  amount  of  depletion.  In 

this  term,  bQ  is  the  fixed  cost  of  depletion  and  j^l  -  F(R)^j  is  the  probability 

of  incurring  the  depletion,  since  F(R)  represents  the  probability  of  issuing 

R  items  or  less  during  the  pipeline  period  when  routine  replenishment  is  occurring, 

and |  1  -  F(R)  |  represents  the  probability  of  demands  being  greater  than  R  during 

that  period.  The  term  (b-j  +  b  o)  J  (x  -  R)  dF(x)  represents  the  expected  costs 

R 

of  depletion  which  vary  either  with  the  quantity  alone  or  with  the  quantity, 
and  duration  of  depletion.  In  this  term,  b-j  represents  the  costs  which  vary 
only  with  uhe  quantity  of  depletion;  b9  p  represents  the  costs  which  vary  with 
both  the  quantity  and  the  duration  of  depletion  multiplied  by  the  length  of 
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%^KS) jJK'V  TEE 
*  *  1  -Lv/  •  .  *'?n 


represents 


'outine  pipe 


Amount  on 


meriataMfOT'cs 
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depletion;  (x  -  R)  is  the  quantity  of  depletion  for  a  demand  x;  dF(x)  is  the 
probability  of  the  demand  x  ;  and  the  integral  from  R  to  oo  represents  the 
summation  of  the  products  of  all  possible  depletion  sizes  and  their  probabil¬ 
ities.  dote  that  the  two  terms  representing  depletion  costs  are  valid  only 
under  assumption  B,  in  which  it  is  assumed  that  the  requisition  is  initiated 
when  the  balance  on  hand  is  exactly  equal  to  the  reorder  level,  R. 

Nov/  we  shall  consider  a  number  of  consecutive  requisitioning  periods, 

0|  (i  a  1  to  n),  But  first,  in  Figure  (?) ,  let  us  join  with  a  straight  line 
the  stock  level  C  at  the  beginning  of  the  requisitioning  period  and  the 
stock  level  D  at  the  end  of  the  requisitioning  period.  This  line  might  be 
considered  as  representing  a  kind  of  "average”  demand  during  the  period  Q< 

If  we  construct  such  lines  for  the  consecutive  requisitioning  periods 
©j,  (i  “  .1  to  n),  we  obtain  the  following  picture: 


From  Figure  (3),  we  can  get  a  clearer  picture  of  our  problem.  If  R  is 
high,  then  we  incur  more  holding  costs.  If  R  is  low,  the  possibility  of  in¬ 
curring  the  iepletion  cost  is  increased*  If  A  is  large,  our  holding  costs 
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are  again  increased.  If  A  is  low,  we  requisition  more  often  and  increase 

the  costs  of  requisitioning.  Our  problem,  then,  is  to  find  the  value  of  R 

1/ 

and  A  which  minimise  the  expected  cost  of  supply  per  unit  of  time.-  To  find 
the  expected  cost  of  supply  per  unit  time  we  must  average  the  cost  of  supply 
over  many  requisitioning  periods  to  allow  for  the  variance  of  one  requisition¬ 
ing  period  from  another.  For  this  reason,  we  cannot  merely  minimize  the  cost 
expressed  by  equation  (l),  since  this  is  only  the  cost  of  one  requisitioning 
period. 

Let  T  represent  a  period  of  operation  of  the  supply  activity,  so  that 

n 

T  a  0.,  +  9_  +  ...  +  9  »  £  0.,  where  9.  is  the  time  period  between  the 

12  n  i-i  1 


i-th  and  the  i+l-th  receipt  of  materiel.  Then  the  total  cost  over  T  is 
n 

given  by  £  L(0. ,R,A),  where  L(9- ,R,A)  is  the  cost  over  the  time  period 
i-1  1 

9j_  as  given  by  equation  (l).  The  cost  per  unit  of  time  is  then: 

n 


cost  per  unit  time  a  X  a 


-  .  n  L(9. ,R,A) 

£  L(9i,R,A)  £  ' - i - 

i°l  _  i“l  n 


n 


0 


E  Gi 

where  we  have  divided  numerator  and  denominator  by  n  and  let  9  »  -  > 

n 

which  is  the  average  requisition  period.  Substituting  from  equation  (l),  we 
now  obtain: 


(2) 


Xs  (d1  +  d^y )  +  1  jk  +  (p1  +  pp  A)  +  (tl  +  t?A)  +  bQ  fl 


F(R)' 


1 J  We  might  remark  here  that  v/e  are  interested  in  minimizing  the  averaged  ex¬ 
pected  cost  per  unit  of  time  rather  than  per  requisitioning  period.  It  can 
be  shown  that  the  value  of  A  which  yields  least  cost  per  requisitioning 
period  is  zero.  This  implies  the  average  requisitioning  period  is  also  zero 
which  in  turn  implies  all  requisitioning  periods  are  zero.  Thus,  we  do  not 
obtain  a  useful  solution  to  our  problem. 
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t  (bl  :  b2P)  Ax  -  R)  dF(x), 

e 

where  we  have  assumed  y  and  A  to  be  constant  relative  to  the  summation* 

• 

In  equation  (?),  y  becomes  the  average  amount  of  stock  on  hand  per 
unit  time,  obtained  by  including  experience  throughout  T.  The  next  thing 
we  must  do  is  express  y  in  terms  of  R,  A,  and  known  quantities.  To  do 
this,  we  introduce  the  notion  of  an  average  requisitioning  period.  The 
length  of  this  period  is  just  9,  of  course.  The  stock  on  hand  at  the  be¬ 
ginning  of  this  average  period  is  R  +  A-  x,  where  x  is  the  average  of  the 
demands  during  the  n  routine  pipeline  time  intervals,  which  immediately 
precede  the  receipt  of  the  materiel;  Similarly,  the  stock  on  hand  at  the 
end  of  the  average  period  is  R  -  x. 

Now  let  the  number  of  requisitioning  periods  become  infinite.  Then 
°o 

x  — =  j  x  dF(x),  which  is  the  average  expected  demand  per  routine 
o 

pipeline  time.  Also,  the  balance  on  hand  during  the  average  requisitioning 
period  will  appear  as  follows: 


Figure  (/+) 

— * 

Amount 

\  '"fL 

r+a-s  ^  J 

I 

L 

i 

on 

Hand 

^~L_ 

Al _  . 

"NS. 

R 

■c - 

-*s. 

R- 

“  — 

• 

0 

Tims 

Notice  that  the  step-function  character  of  Figure  (A)  is  due  to  the  fact  that 
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we  issue  one  item  at  a  time,  and  cannot  issue  fractions  of  items.  If  the 
amount  on  hand  were  a  linear  function,  as  represented  by  the  diagonal  dotted 

line  in  Figure  (4),  the  average  amount  on  hand  would  just  be  R  -  £  + 

, ,  a  term  of  1  to  allow  for  the  step-function 

However ,  we  must  add  a  correctio  ^  p 

effect.  Thus,  we  obtain 
(3)  y  =  R  -  €  *  ^  J 

The  value  for  y  in  equation  (3)  validates  our  prior  assumption  that 
y  be  constant  relative  to  the  suction  in  equation  (?).  It  might  also  be 
noted  that  we  are  charging  the  warehousing  costs  for  only  an  amount  y  whereas 
warehouse  space  is  needed  for  the  maximum  amount  stocked.  We  feel  our  assump¬ 
tion  justified  on  the  basis  that  all  items  in  a  warehouse  will  not  be  stocked 
in  their  maximum  amounts  at  the  same  time,  but  will  indeed  average  out  as 
assumed,  with  some  items  requiring  space  for  more  than  the  amount  assumed  and 

others  less  in  any  given  period  of  time. 

Next  we  will  establish  a  relationship  between  A  and  9.  If  we  divide  T 

into  n  equal  intervals,  each  oflength  9,  we  can  write 

0 

Z  x, , 

,  ,  ik 
x .  3  k“l 

1  — — - 


as  the  arithmetic  mean  demand  for  the  i-th  period.  Therefore, 

9  .  t 

a  ~  a  Ex.,  and  summing  over  1,  we  geu 
w  i  k-1  ik 

1/  A  -re  rigorous  development  of  equation 

ffirtSTwith \£  Ihen  the  average  stock  on  hand  during  the  k-th  day 
of  0.  is  given  by 

J  k  /,  ^  ^  \ 


FiK  ■  R  -  E.XiJ  *  ^  +  5 
»  )  ^ 


(1  <  k  <  p) 


where  we 


R  *  4  .*  t^ij  *  f ik  *  I  >  (p  <  k  <  p  9i) 

■si.  mu™  •<  renresent  a  day's  demand  rather  than  the 


are  temporarily  letting  x  represen 
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n  _  n  Q  n  ®i 

0  £  x4  3  Z  £  x  “  £  £  x  =  n  A, 

i“l  x  i-1  k-1  ik  i=l  j»l 

which  is  the  total  demand  during  T.  Employing  the  extreme  left  and  right 
members, 

n  _ 

£  x 

(4)  A  "  ©  ial  1  •  3  0  x. 
n 

Estimating  x  by  £  ,  we  substitute  the  values  for  y  and  0  from 
equations  (3)  and  (4),  respectively,  into  equation  (?)  and  get 


X(A,R)  «  d1  +  d2(R 


f  + 


e(k  +  P  -t1  +  b0> 

A  +  1)  + _ _ _ _ _ 

2  A 


+  £(p?  +  t.?) 


F(R)  + 


A 


(x  -  R)  dF(x). 


demand  per  pipeline  time  as  previously  defined.  To  show  the  derivation  of 
these  equations  more  clearly,  we  might  isolate  the  k-th  day,  which  would 
appear  as  follows: 


R  - 


Notice,  in  particular,  the  correction  factor  l/?  in  the  diagram.  In  general, 
if  we  issue  m  items  per  issue  and  n  issues  per  day,  equally  spaced  throughout 
the  day,  then  the  amount  on  hand  during  the  time  period  between  the  i-th  and  the 
i+l-th  issues  is  just  mi  (i  *  n,  n  -  1,  ...,  1)  plus  the  amount  on  hand  at  the 
end  of  the  period.  The  average  amount  on  hand  during  the  day  is  then 


7 
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This  equation  can  be  rearranged  to  get 


(5)  Xa, R)  "  A  ♦ 


e(B  +  bo)  d2  __  fb 


+  .  +  d2R  -  F(R)  +  fC  J  (x  -  R)  dF(x), 


A  R 


where  A  ®  d1  +  £  (p^  +  t^)  +  dy(l  -  €  ) 


B  e  k  +  +  t^ 


C  *  +  b^p. 

The  function  X  is  minimized  by  equating  to  zero  the  partial  derivatives 
of  X  with  respect  to  A  and  R  and  then  solving  these  equations  simultan¬ 
eously  to  provide  optimal  values  for  A  and  R.  Therefore,  we  obtain 
6  b  n°° 

(6)  ”  d  -  r(R)  +fc  a  ,/  (x  -  R)  dF(x)  «  0 

OR  A  T  £R  °R 


E  mi  »  mn(n+l)  ®  mn  +  m  (plus  amount  on  hand  at  the  end  of  the  period). 

ial  n  2n  2  2 

Hence,  the  correction  factor  is  m.  Under  assumption  A,  of  course,  ra  a  1  and 

2 

the  correction  factor  is  1.  From  the  equations  (3a),  we  get  the  average 

2 

amount  on  hand  during  9^  to  be: 


(e  -  l)  Pyi 

y .  3  R  +  l  + _ a  -  ^ 

1  2  9.  j-1 


pei  (pe,  -  j  +  l) 


?  -  J  +  1)  x  P®i  X 

^ _ W  +  1  £  _ 

p0^  2  k"l  p0^ 


(3b)  y, 


P®!  (j  -  1)  x 

R  +  i  -  _A  +  E  2 

2  9.  j3!  p0. 

i 


where  we  have  set  A  “  E  x.  ..  However,  by  so  expressing  a*  we  infer  that 

j-1  1J 

the  balance  on  hand  goes  negative  in  response  to  demands  after  total  stock 
depletion.  In  actuality,  balances  do  not  go  negative;  but  the  correction  term 
to  be  added  to  the  value  for  y  which  we  obtain  by  allowing  negative  balances, 
can  be  shown  to  be: 


d2  f(B  +  bn)  -  f bnF(R)  +  fC /  (x  -  R)  dF(x) 

(7)  9a  " _  - _  _ _ _R _ _ 53  0 

aA  2  " 
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We  may  solve  explicitly  for  A  in  equations  (6)  and  (7)  and  get 

(8)  A  =  i_U  f (R)  -  C  A  (x  -  R)  dF(x) 

dJ  “  »  J  R  3 

>  .  v 


HR 


(9)  A< 


roo  7 

26  f  B  +  bn  [l  -  F(R)]  +  C  /  (x  -  R)  dF(x)  f  . 
do  1  "  R  J 


Now  if  A  is  eliminated  from  equations  (8)  and  (9),  the  optimal  reorder 
level,  R#,  may  be  determined  from  the  equation  which  results  and  substituted 
back  in  equation  (8)  or  (9)  to  find  ,  the  optimal  amount  of  reorder. 
Whether  or  not  the  values  so  obtained  yield  the  actual  minimum  supply  cost 
depends  upon  an  investigation  of  the  second  partial  derivate3.  We  will  omit 
this  investigation:  for  all  reasonable  demand  functions,  a  valid  minimum  does 
exist. 


£  f 

A  R 


[x  -  R)  (x  - 

4  x 


dF(x). 


This  term,  upon  inspection,  is  so  small  that  it  can  be  comfortably  neglected. 
Returning  to  equation  (3b),  we  take  the  weighted  average  over  n  requisitioning 
periods  and  obtain  ,  _ _  > ,  a  ~ 

j  9  v.  n  (j  -  1)  x  yti 

‘  A  1-A+LE  2  V-4'* 


(3c)  y  ■  ial 


iral  j*l 


By  re-subdividing  our  interval  T  into  n  equal  intervals,  each  of  length  6, 

we  can  rearrange  terms  in  the  last  member  of  the  above  equation  to  get 


n  P0!  (j  -  1)  x  n  p9 
v  I  2  ^  "PS 

i-l  3-1  „„ff  i”1  k“1 


1)  x  n  pO  (k  -  i)xik  pQ  (k  -  l)  i 

2  J  ? _  mV  ?  Am 


n  x., 

rt  lk 


2  i»l  n 


p9  (k  -  1)  >\ 


But  x  e  for  each  k  as  n  — ^  oo, 
P 


Substituting,  we  got 


/ 
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The  calculation  of  R#  and  A*  from  equations  (8)  and  (9)  is  a 
relatively  easy  matter,  particularly  on  the  electronic  calculators.  In  fact, 
for  some  simple  demand  functions,  explicit  answers  for  R*  and  A#  can  be 
obtained.  Examples  of  such  functions  are  discussed  in  a  later  section. 

The  case  of  recoverable-type  items  at  the  base  can  be  included  in  the 
above  results  by  just  considering  the  demand  function,  f(x),  to  represent 
the  net  loss  due  to  items  condemned  and  items  beyond  base  repair.  This  is 
due  to  the  fact  that  if  an  item  is  base-reparable,  it  can  readily  be  con¬ 
verted  to  a  serviceable  item  -  more  readily,  in  fact,  than  obtaining  the 
item  from  any  other  source.  Hence,  a  base-reparable  item  can  be  treated  as 
s  ervieeable  from  a  stock  policy  point  of  view. 

In  the  derivation  of  our  base  results  as  expressed  by  equations  (8)  and 
(9),  we  used  the  idea  of  the  number  of  requisitioning  periods  becoming  infinite. 
One  might  question  whether  the  results  are  tenable  based  upon  such  processes, 
inasmuch  as  it  is  certain  the  supply  system  will  not  operate  forever.  Actually, 

n  P°i  (j  -  l)x  p©  (k  -  l)^ 

^  2  2'  1J  »  £  2  "  » 

i“l  j-1  _  k®l 

np0  p'  9  2p  0 

*  fj 

2 

Substituting  back  into  equation  (3c),  we  get 

(3d)  y  “  R  +  1  -  A_  +  £  Q  . 

2  -  2 

0 

Subsequently,  we  show  that  A  *  £  9  so  that  equation  (3d)  becomes 

y  “  R  -  £  +  411  . 

2 
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however,  the  method  is  validated  by  an  implicit  assumption,  namely  that’  our 
demand  probability  function  is  quite-  independent  of  time.  Thus,  even  though 
reason  might  assess  a  zero  probability  of  issuing  n  items  a*  thousand  years 
from  now,  this  probability  is  not  reflected  in  our  demand  function.  Indeed, 
our  demand  function  assesses  the  probability  of  demand  to  be  the  same  for  • 
all  time.  But  since  our  problem  was  to  establish  some  A  and  R  to  use  in 
the  stockage  policy,  the  use  of  A*  and  R*  as  calculated  above  is  the  best 
we  can  do  if  the  given  demand  function,  f(x),  is  the  best  estimate  of  demand 
probability  that  we  can  obtain.  If  the  demand  probability  function  can  be 
expressed  in  terms  of  time,  then  this  becomes  another  problem  for  which  we 
have  no  general  solution.  On  the  other  hand,  we  can  approximate  a  solution 
by  recalculating  R*  and  A*  periodically  with  a  demand  function  adjusted 
to  reflect  the  trend  established  by  past  consumption  and  other  factors.  The 
derivation  of  the  demand  probability  functions,  however,  is  beyond  the  scope 
of  this  paper. 

3-  ROUND IMG  TO  INTEGRAL  VALUES 

For  practical  application  of  the  optimal  reorder  level,  Rw  ,  and  reorder 
amount,  A#  ,  we  must  examine  the  problem  of  rounding  them  to  integral  values. 
If  we  construct  our  cost  surface,  X(A,R),  in  the  vicinity  of  (as,R*),  we  might 
obtain  the  situation  shown  in  Figure  (5). 

Figure  (.5) 


/ 
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In  Figure  (5),  R*  and  A»  represent  the  largest  integers  contained  in 
R*  and  A#  respectively.  The  values,  Rr,  R*  +  1,  A1,  A,+  l,  determine 

the  rectangle  ABCD  on  the  cost  surface.  At  first  glance,  it  might  seem 
justified  to  round  A#  and  R#  to  those  integral  values  of  A  and  R 
which  yield  the  least  cost  among  the  four  values,  l(RT,  A*)»  XR^A'+l), 
XR,+1,A‘),  and  XR,+1> A'+l)  which  occur  at  the  points  A,  B,  C,  and  D.  On 
the  other  hand,  it  is  entirely  conceivable  that  the  cost  surface,  X(A,R), 
in  the  vicinity  of  (Att,R#)  might  appear  as  follows: 

Figure  (6) 


From  Figure  (6),  it  is  clear  that  it  is  quite  possible  for  the  least  cost  X, 
for  integral  values  of  A  and  R,  to  occur  not  at  the  corners  of  the 
rectangle  ABCD  containing  >(  A*,R#),  but  at  the  points  E  or  F  or,  in  some 
cases,  at  even  more  remote  points  on  the  surface. 

In  general,  then  to  find  the  integral  values  for  A  and  R  which  afford 
the  least  cost  of  supply,  we  must  proceed  as  follows:  First  we  find  that 
value,  A1  ,  which  minimizes  the  value  )(R  *  Rf ,A).  We  do  this  by  substituting 
R<  for  R  in  equation  (7),  the  partial  derivative  of  X  with  respect  to  A , 
and  solve  for  A  ,  the  solution  being  A^#  .  Wc  then  repeat  this  process  by 
substituting  R»  +  1  for  R  in  equation  (7)  to  find  A2>  that  value  of  A  which 


minimizes  X(R  =  R?+1,A).  Continuing  this  procedure,  we  substitute  A1  for 
A  in  equation  (8),  the  partial  derivative  of  X  with  respect  to  R,  and 
solve  for  R^,  the  value  of  R  which  minimizes  X(R>  A=  A*)*  Similarly, 
we  find  R^,  the  value  of  R  which  minimizes  X(R,  A  “  AT  +  l).  We  have 
now  found  the  points  (R*,  A^),  (Rf  +  1,  A^),  (R^,  AT),  and  (R?,  Af  +  l)  on 
the  X  -  surface,  which  represent  the  minimum  values  for  X  along  the  ex¬ 
tended  "sides”  of  the  rectangle  ABCD  in  Figure  (5).  If  we  let  the  symbol 
H  represent  the  largest  integer  contained  in  x,  we  can  now  set  down  the 
following  eight  values  for  X  J 


(cp 

Xn<,  gy  ) 

(c5) 

X(  Dp],  A') 

(C?) 

mt',  [>.,>1) 

<V 

xDpli,  a') 

(CI3) 

Xr'+i,  Dy )' 

(°?) 

iOjl  A'+l) 

oy 

Xrui,  [a?1D 

(V 

xCRpll.AHl) 

Of  the  eight  values  of  X  so  obtained,  the  least  one  provides  the  integral 
values  of  A  and  R  to  be  used.  The  conditions  on  the  nature  of  the  surface, 
X(A,R),  which  validate  this  process  are  assumed  to  exist  and  do  exist  for 
any  practical  case. 

It  might  be  remarked  that  the  situation  depicted  in  Figure  (6)  can  occur 
in  practice  only  for  low  values  of  A#  and  Rtf,  For  large  values  of  A# 
and  Rtf,  the  cost  surface  becomes  so  flat  that  one  can  automatically  round 
to  the  nearest  integer.  For  a  region  of  the  surface  between  these  extremes, 
only  the  four  values,  X[Rf,Af),  Rf+l,Af),  X(Rf,At+l),  and  x(Rt+l>Af+l) 
need  be  computed  and  compared. 


4.  DETERMINING  WHEN  TO  STOCK  AT  THE  BASE 


Since  our  cost  function,  X(A,R)  in  equation  (5),  expresses  the  cost 
of  holding  an  item  of  supply  at  the  base,  it  is  invalid  for  A  *  0,  R  *  0. 
This  situation  would  be  interpreted  as  not  stocking  the  item  at  the  base  at 
all,  but  stocking  the  item  at  the  depot  instead.  The  average  cost  of  supply 
per  unit  time,  in  this  case,  becomes 

(xo)  l\  -  £  (b  +  b:  +  b?p) 

which  is  nothing  more  than  the  depletion  penalty  multiplied  by  the  average 
probability  of  incurring  it  per  unit  of  time.  A  superposition  of  the  two 
cost  surfaces  represented  by  equations  (5)  and  (10 )  might  appear  as  follows: 

.  i 

Figure  (7) 


In  Figure  (7),  the  hatched  area  of  the  A-R  plane  represents  those  values 
of  A  and  R  for  which  the  X  -  surface  lies  below  the  ^  -  plane.  From 
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equation  (5),  we  also  notice  the  cost  X  approaches  infinity  as  A  approaches 
zero,  which  renders  the  X-  surface  invalid  at  (0,0).  Of  course  this  situation 
is  reasonable;  if  A  approaches  zero  so  does  0  and  the  fixed  costs  for  each 
0^  occur  more  and  more  often  until  in  the  limit  they  occur  infinitely  often. 

In  Figure  (7),  it  is  clear  that,  with  the  scale  values  chosen,  the  cost  ■ 

A  on  the  X  -  surface  (which  is  yielded  by  A  =  1,  R  **  0)  would  be  chosen 
for  any  optimal  A#  and  R*  in  the  cross-hatched  area.  A  very  simple  method, 
however,  can  be  established  to  decide  whether  to  stock  at  the  base  or  not. 

This  method  is  to  compute  the  optimal,  rounded  A  #  and  R#  as  outlined  in 
sections  2  and  3,  and  then  compare  the  resulting  cost  A  with  the  cost,  t[  ?  . 
of  not  stocking  at  the  base.  We  notice,  in  this  connection,  that  the  cost 
is  easy  to  calculate,  being  nothing  more  than  the  cost  of  depletion  multiplied 
by  i  . 

From  our  expression  for  X  in  equation  (5),  it  is  fairly  clear  that  the 
possibilit3>-  of  not  stocking  at  the  base  can  occur  only  for  high-cost,  low  de¬ 
mand  items  or  for  low  demand  items  with  low  depletion  penalties.  For  low-cost, 
low  demand  items,  the  constant  costs  of  requisitioning,  in  addition  to  any 
relatively  high  depletion  costs,  will  cause  the  items  to  be  stocked  at  the 
base.  In  fact,  a  high  depletion  cost  for  low-cost  items  will  cause  large 
stocks  at  the  base  even  for  very  low  demand  rates. 

5.  OPTIMAL  POLICY  FOR  UNKNOWN  MEAN  DEMAND 
If,  in  the  results  of  section  2,  the  form  of  the  density  function  f(x) 
is  known  but  the  mean  demand  €  can  be  expressed  only  as  a  probability 
function,  g(t),  we  can  extend  our  theory  as  follows:  First  we  calculate  our 
optimal  cost  as  a  function  of  €  : 

•(31)  \(£)  -  XfA*,R*,£)- 

Then,  for  each  optimal  A*  and  K-  we  multiply  the  cost  x(0  by  the  prob- 
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ability  of  incurring  that  cost  and  sum  over  all  such  probabilities.  In  this 
manner  we  obtain 
roo 

(12)  r(6)  «  J  g(t)  \A*(£),R*(6),t)  dtp 

o 

in  which  we  changed  the  variable  £  to  t  for  the  integration  and  remembered 
that  A*  and  R#  are  themselves  functions  of  6,  That  value  for  6,  then, 
which  gives  least  cost  is  calculated  by  the  usual  procedure  of  setting  the 
first  derivative  equal  to  zero.  Therefore,  optimal  A*  and  R#  are  given 


A#  »  A«(g#) 

(13) 

R#  »  R#(£#), 

where  6#  is  a  solution  of  the  equation 


dV(fc)  B  0  *»  _d_  f  g(t)  x(A»(6),R^(€),t)  dt. 
d  €  d6  o 


The  rounding  procedure  for  A*  and  R^,  in  this  case,  becomes  rather 

cumbersome.  In  theory,  the  procedure  would  proceed  as  follows:  Let  R*  -Of! 

(the  largest  integer  contained  in  R#),  as  previously  defined.  Proceeding  as  in 

section  3,  we  obtain  A^*(£^)  as  a  solution  of  the  equation  0X(R* ,  A, 6)  =»  0. 

5A 

We  then  follow  the  procedure  in  the  first  part  of  this  section  to  obtain 


Y(fx)  ”  J  g(t)  MR',4^  (fp.O  dt, 

o 

and  then  derive  as  a  solution  of  the  equation 
manner  we  get  the  costs 


dY(f  J 


0.  In  this 


(C,)  X(R',A1)  am!  (c2) 
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where  A^  -  [~A^  (6^)]] 

?  OT  (6*0  +.1. 

By  a  similar  process,  we  can  derive  the  costs 
(C3)  X(R»+1  »/^)  and  (C/+)  XRV+1,  A^), 
where  A^  =  [a^  JO 

^4  *  \LL2^p~]  +  1 

A  t  =  the  solution  of  d.X(Rf+l,At6 )  =  0. 

dA 

Through  this  method  we  can  also  get 

(C  )  XO^.A') 

(c6)  x(R?,a') 

(C7)  X(R  ,a*+i) 

(cg)  x(R^,a'+i) 

wherein  the  roles  of  A  and  R  have  been  reversed*  The  eight  cost  values 
so  obtained  can  then  be  compared  to  find  those  integral  values  of  A  and  R 
which  yield  least  cost. 

In  practice,  most  of  this  procedure  would  be  done  by  numerical  methods, 
particularly  if  explicit  values-  for  R**  and  A*  could  not  be  obtained  from 
equations  (8)  and  (9).  If  this  were  the  case,  even  R*  and  A#  for  known  6 
would  probably  be  found  by  numerical  methods. 

One  such  numerical  method,  for  6  given  by  g(t),  might  be  as  follows: 

(a)  Divide  the  interval  £a,bl  on  the  t-axis  of  the  function  g(t)  into  n 
equal  sub  intervals,  where  the  points  a  and  b  might  represent  confidence 
limits  at  some  level.  Let  6  .  (i  a  1  to  n)  represent  the  midpoints  of  the 
subintervals  so  obtained,  with  g(6^  )  denoting  the  associated  probability. 

(Thus  we  approximate  the  continuous  function,  g(t),  by  a  density  function, 
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where 


g(6^)  =  J g(t)  dt,  and  ^  is  the  i-th  subinterval). 

(b)  For  each  compute  optimal  rounded  A  and  R  as  in  section  2  and  J, 
Let  these  optimal  values  be  denoted  and  A?. 

(c)  Compute  the  matrix  X  ^  Form  the  matrix  product 


£  \  ±  g(£'J  s  Yj  •  (i  8  1  tp  n) 
j=l  J  J 

(d)  Let  Yyf.  represent  the  least  of  the  values,  Y^*  Then  A^,  and  R#-, 
provide  the  optimal  values  for  A  and  R. 

6.  EXAMPLES 

To  illustrate  previous  results,  we  will  now  consider  three  examples  of 
particular  demand  probability  functions.  In  the  first  example,  we  are  able 
to  get  explicit  answers  for  both  known  and  unknown  mean  demands.  In'  the 
second,  we  can  get  explicit  answers  only  for  the  case  of  known  mean  demand, 
while  in  the  last  example,  we  can  go  scarcely • further  than  we  have  already. 
Example  1.  Let  the  demand  probability  function  be  of  the  form: 

f(x)  ■  a  x  =  0 
(14)  =  P  x  *»  1 

=?  0  all  other  x. 

This  demand  function  is  very  useful  to  express  the  case  of  very  low  demands 

providing  we  take  (3  quite  low,  say  less  than  0.1. 

Our  condition  that  j  f(x)dx  =  1  implies  that  a  +  p  «  1,  or  (3  “  1  -  a. 

o 

Also  we  see  that  oo 

£  “  P  x  dF(x)  »  p 
^ o 

Since  the  last  two  terms  of  the  cost  function,  equation  (5),  remain  the 
same  for  all  R  >  1,  it  is  clear  that  minimum  cost  can  occur  only  for  R  ■  0,  or 


■trf-tffir.il 


A 
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r  *  1  only  integral  values  for  R  are  allowed.  By  direct  substitution 
in  equation  (5)#  we  get 

dP 

X(A,R«0)  -  A  +  _  A  +  p(?  ,tpHl 
2  A 

(15) 

X(A,R=l)  «  A  +  d  (  A  +  l)  +  j3T 
^  2  A 

where  T  -  k  +  P-j_  +  ^  and  TT  ■  bQ  +  C  =>  bQ  +  bj  +  b2p  =  total  depletion  cost. 
Equation  (9)  yields  the  values 


2j  (T  +  {371} 
d0 


(using  R  B  0) 


A* 

1 


(using  R  =  l) 


It  is  clear  that  the  equation  X(A*fR-0)  -  X(A|,R»l)  .yields  a  p»  such 
that  for  p  <C  p  * ,  R  “  0  with  A^  provide  least  cost,  and  for  p  2*  p  ,  R  1 

with  A*  provide  least  cost.  A  calculation  shows,  in  fact,  that  pf  must 
satisfy  the  equation 


Q2p»?7T-  dp]  2  =>  8  d2p»T. 

For  rounding,  we  compare  the  four  values 

xCa$  Ra°) 

X([A^1,  R-0) 

xOfl  R=,i) 

X([Aj>l,  R“l) 

and  choose  the  integral  A#  and  R«  which  provide  the  least  value.  The  least 
value  so  obtained,  X(A«,R«)  car.  then  be  compared  against  the  cost  f|  - 
p(b  +  b  +  b  p)  to  decide  whether  or  not  to  stock  at  the  base,  *v'e  would  not 
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stock  if  7\^  <  A(Atf,R*). 

Let  us  now  consider  the  case  where  p  ,  the  arithmetic  mean  of  f(x)#  is 
given  by  a  probability  function,  g(t).  Equation  (12),  together  with 
equations  (15)  give 

,<fl  "  *  *  ?  W 1  *  (H  -  .) 


A*(p) 

r(p)  =  A  +  d2( _  + 1)  + 


.jgT 


(R  -  l) 


where 


-  z1 00 
g  83  /  t  g(t)dt 
o 


P  o 


oo 


t  g(t)dt. 

The  derivatives  with  respect  to  {3  yield 


-  SLltl f 


A#  2 
0 


gT 

2  Ai2 


3A» 

■dp 


^1 

ap 


from  which,  with  squationa  (16),  we  find 


(R“0) 


(R  -  1) 


=*  -  T  +  Vr  -»•  4^  (gT  +.U.1P)  (R  *  0) 


7f 


}i  -  s 


wherein  we  have  ignored  negative  values  for  p#  .  We  then  calculate  the  four 
values, 
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(cl) 

X(&0,R=0) 

(°2) 

X(A  +1,R»1) 

(C3) 

^  ,R-1) 

(V 

where  *  Q^q(Pq)D 

ax  -  k  (pja 


the  least  of  which  provides  optimal  integral  values  for  A  and  Pw.  The  cost 

so  obtained  should  be  compared  against  the  cost  of  not  stocking  at  the  base, 

which  is  given  by  n  =  g  (b  +  b  +  bn)  *  g  7T 

0  I  - 

Example  2.  Let  the  demand  probability  function  be  of  the  form 
(18)  f(x)  ■  ae 

This  function  is  useful  for  somewhat  higher  demand  rates  than  in  the  previous 

example,  and  is  also  useful  in  illustrating  the  procedure  for  the  case  of  a 

continuous  demand  function. 

00 

Our  condition  that J  f(x)dx  a  1  implies  that  a  °p.  Also  we  calculate 


i  -  f*  x  f(> 


f(x)dx  to  be  1. 


a 


Substituting  f(x)  a  ae“aA  and  F(x)  *  j  f (t)dt  a  1  -  e  ^  into 


equations  (8)  and  (9),  and  eliminating  A  ,  we  obtain 

2d. 


( 


-aR* 


(K  *  C  )  1  ‘ 


r  -clR#  /. 
e 


a 


(bQ 


Let 


e~cR*  (b  +  C).  Then 
0  a 

2d„  2d 


z  - 


a 


f  «  -  2  B  -  0, 

a 


and 


r 


1  + . 


1  +  2a  B 

u 
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Since  z  is  positive  and  all  terms  under  the  radical  are  positive,  we  take 
the  plus  sign  for  the  radical  and  get 


1-  .  V  *  ' 


+  Pa  B 
d~ 


having  replaced  our  value  for  z.  When  substituted  back  into  equation  (8), 
we  get 


A*  - 


e  (bQa  +  C) 


1  +  q/l  +  2aB 

*  do  . 


Also,  the  optimal  average  requisitioning  period  is  given  by 
9*  =>  Af  ®  aA,Uf  “  1  +*|/  1  +  2a  B. 

e  *  d2 

We  can  summarize  our  results  by  the  equations 


9*  *»  1 


+  V1  + 


Pa  B 

dr. 


A-  =  0* 


P*  =*  1  lr 


V  +  c 


d  0» 
2 


Integral  values  for  A#  and  R  #  can  now  be  found  by  the  method  of 

section  j.  For  the  situation  where  the  mean,  1,  can  be  expressed  only  by 

a 

a  probability  function,  g(t),  an  explicit  answer  for  optimal  a  and  R 
cannot  be  obtained;  in  practice,  a  numerical  method  would  be  used  as  suggested 


in  section  5» 


Example  3*  Let  the  demand  probability  function  be  of  the  form 


(20)  f(x)  *  e~mmx  (Poisson  distribution). 

x  J 

The  Poisson  function  is  useful  to  express  most  demands  at  bases,  provided 
it  is  used  on  an  item  of  supply  basis  as  is  done  in  this  paper.  This ' function, 
in  particular,  is  much  more  realistic  than  either  of  the  other  two  examples 
for  higher  demand  rates. 

00 

We  see  that  this  function  satisfies  our  condition  that  J  f(x)  dx  *  1 

..po  0 

and  we  also  calculate  6  j  x  f(x)  dx  to  be  just  m.  For  the  Poisson 

"0 

function,  equations  (8)  and  (9)  become 


(21)  A*  -  _m 

d 


2 


|  bQ  f(R»)  +  C  f  1  -  F(R*-1)  | 


(22)  A*  -  2m  j  B  ♦  \  +  C  (m  -  R 


R»)  “l  -  F(R*)“  +  mC  f(R#)  j 


Eliminating  A  ,  we  get 


(23) 


h(R»)  -  m|b0f(R»)  +  C  f  1  -  F(R«-1)  j  2  -  +  p>0  +  C(m  -  R»)J 


1  -  F(R«) 


+  mC 


0. 


This  equation  can  be  solved  for  by  numerical  methods  and  R#  substituted 
back  into  equation  (21)  for  Atf  .  The  solution  of  equation  (23)  provides- two 
integral  values  for  R*,  one  of  which  yields  the  least  positive  h(Rtt)  and  the 
other  the  greatest  negative  h(R*).  These  two  values,  R|  and  R#,  when  sub¬ 
stituted  into  equation  (21),  yield  A*  and  A# .  These  values  for  A  are  then 
rounded  to  QdjJ#  D^J+1,  l_A!p»  G'llr  Four  values  for  X  are  computed  as 


follows : 
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(i)  XD£L  R|) 

(3)  R|) 

(2)  Xp£K  HJ) 

(4).  xO|J+l,  H8.) 

and  those  integral  values  for  R  and  A  which  provide  the  least  of  the  four 
c  osts  are  chosen  for  our  optimal  A  and  R. 

For  m  expressed  as  a  probability  function,  g(t),  the  solution  must 
be  obtained  by  numerical  methods  as  outlined  in  section  5» 

7.  MODIFICATION  OF  THE  BASE  COST  .EQUATION 
The  base  cost  equation,  as  expressed  in  equation  (5),  may  be  modified 
and  refined  in  several  respects.  One  such  refinement  has  already  been 
suggested  in  footnote  3* 

Another  correction  to  the  equation  should  be  made  to  allow  for  the 
possibility  of  a  routine  stock  replenishment  arriving  while  a  priority  requi¬ 
sition  is  being  processed.  Thus,  if  the  priority  pipeline  time  is  2  days, 
the  cost  equation  as  expressed  charges  a  7  day  depletion  penalty  for  a  pri- 
ority  requisition  submitted  1  day  before  the  receipt  of  the  routine  replenish¬ 
ment.  Of  course,  only  1  day*s  depletion  penalty  should  actually  be  charged. 
The  correction  factor  to  allow  for  this  overcharge  is 

b  pp  aOo 

'  J  (x  -  R)  dF"(x),  where  FM  is  a  demand  probability  function 
2  R 

expressed  in  terms  of  p  -  pas  the  unit  of  time.  If  f(x)  is  a  Poisson  dis¬ 
tribution  expressed  as 

f(m,x)  a  e~nimx  , 

xl 

then  Fn(x)  is  given  by 

F"(x)  «  /  f(  (£ip)m,t)  dt. 

'  o  p 

This  correction  factor,  unlike  that  mentioned  in  footnote  3,  may  be  of  appre¬ 
ciable  magnitude. 


P-647 

-34- 


Another  refinement  in  the  base  cost  equation  may  be  made  by  considering 
an  addition  assumption: 

F.  A  priority  request  is  not  submitted  for  any  depletion  occurring  in  the 
last  p  days  of  the  routine  pipeline  time. 

Under  this  assumption,  the  base  cost  equation  is  further  amended  by  sub-  • 
stituting  the  function  Fu  for  F  in  the  bQ  and  b^  terms  of  eouation  (5)« 

Another. way  in  which  our  cost  equation  may  be  modified  under  further 
assumptions  is  in  the  manner  of  submitting  priority  reouisitions.  Assumption 
A,  which  causes  demand  to  occur  one  item  at  a  time,  together  with  assumption 
D,  more  or  less  implies  that  a  priority  requisition  be  submitted  for  each  item 
demanded  that  is  not  available  at  the  the  base.  In  this  case,  our  term  b  is 
zero  (unless  we  wish  to  assess,  in  some  way,  the  part  of  the  capita,  costs 
necessary  to  establish  facilities  for  the  special  handling  of  priority  requi¬ 
sitions,  whether  or  not  such  facilities  are  actually  used  in  a  particular  case), 
and  the  term  b^  includes  all  costs  of  paperwork,  handling,  communications, 
packaging,  inspection,  and  transportation  for  the  priority  request.  Each  item, 
of  course,  is  packaged  and  shipped  separately.  However,  it  would  seem  more 
reasonable  to  institute  another  level,  Rp,  less  than  or  equal  to  Rj  such 
that  when  this  level  is  reached,  a  priority  requisition  be  submitted  for  an 
amount  A  .  The  priority  reorder  level,  R^,  may  range  from  R  down  to 
3ome  negative  value.  This  concept,  in  fact,  may  be  extended  to  provide  as 
many  reorder  levels  and  reorder  amounts  as  there -are  modes  of  transportation, 
and  would  in  theory  exactly  specify  the  kind  of  transportation  each  shipment 
to  the  base  would  receive.  However,  to  apply  such  a  concept  would  require 
that  the  costs  and  pipeline  time  for  each  mode  of  transportation  be  constant; 
this  requirement  of  course  is  unrealistic. 

A  more  realistic  approach  toward  the  base  stockage  policy  would  be  to 


< 
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recognize  two  types  of  requisitions  as  represented  by  two  general  modes  of 
communications  and  transportation.  One  type  of  requisition  would  be  the 
routine  requisition  instituted  when  the  amount  on  hand  reaches  the  reorder 
level  R.  The  pipeline  time  for  this  kind  of  requisition  should  be  con¬ 
sidered  as  a  random  variable  to  allow  for  the  vagaries  of  slower  communica¬ 
tions  and  surface  transportation.  The  other  type  of  requisition  would  be 
the  priority  requisition  instituted  when  the  amount  on  hand  reaches  the  re¬ 
order  level,  Rp.  The  pipeline  time  for  this  kind  of  requisition  may  be 
considered  constant  inasmuch  as  electrical  communications  and  air  trans¬ 
portation  yield  much  less  variable  pipeline  times,-  In  this  manner,  we 
derive  two  possible  modifications  of  our  base  cost  equation: 

(1)  Inclusion  of  a  priority  reorder  level,  R  ,  and  a  priority  re- 

P 

order  amount,  A  . 

P 

(2)  Expression  of  the  routine  pipeline  time  as  a  random  variable  in¬ 
stead  of  being  constant. 

The  first  modification  is  being  subjected  to  further  study.  The  second 
modification,  however,  can  be  readily  accomodated  by  the  procedure  of  section 
5,  providing  the  demand  function  is  Poisson  distributed,  A  characteristic 
of  the  Poisson  distribution  is  that  the  probability  of  demand  is  proportional 
to  the  length  of  time.  Therefore,  we  may  express  6  as  g(t)m,  where 
g(t)  is  now  interpreted  a3  the  probability  function  for  the  routine  pipeline 
time,  expressed  in  days,  and  m  is  the  expected  mean  demand  per  day.  This, 
of  course,  phrases  the  problem  exactly  as  in  section  5-  Of  course,  vre  can 
still  allow  for  the  case  of  unknown  mean  demand  in  addition  to  random  routine 
pipeline  time  by  expressing  m  as  rr(s),  so  that  €  is  given  by  g(t)  m(s)  . 
The  procedure  of  section  5  still  applies  directly,  but  with  double  integration 
since  the  variables  s  and  t  are  independent. 


/ 
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8.  APPLICATION  OF  THE  OPTIMAL  INVENTORY  POLICY  TO  A  COMMERCIAL  FIRM 

The  foregoing  formulation  of  an  optimal  inventory  policy  is  applicable 
to  a  profit-making  firm  with  only  minor  modifications.  We  redefine  our  costs 
as  follows: 

(a)  d„  +  d  y  a  the  cost  of  holding  in  inventory  an  amount  y  for  one 

5  4 

unit  of  time,  including  depreciation,  obsolescence, 
warehousing  costs,  interest  costs,  and  insurance. 

(b)  p^  +  p^  +  p^A  »  the  cost  of  buying  an  order  of  size  A, 

(c)  bQ  +  cq  =  the  cost  of  depletion  by  amount  q,  including  all  present 

and  future  losses  of  custom,  discounted  to  the  current 
time  period. 

(d)  M  »  the  average  gross  utility  of  operating  the  inventory  for  one 

unit  of  time, 

(e)  All  other  symbols  retain  their  former  definitions. 

We  may  now  write  the  equation  for  the  net  loss,  L,  during  a  requisi¬ 
tioning  cycle,  0,  as: 

L(R,A,e)  “  (d  ♦  d  )  9  +  k  +  (Pl  +  pA  ♦  p^a2)  +  (tx  *  t?A)  +  bQ  |JL  -  F(R)] 

n  OO 

+  c  J  (x  -  R)dF(x)  -  MO. 

R 

Since  all  costs  and  utilities  in  each  requisitioning  cycle,  0,  are 
discounted  to  the  current  time  period  (c  was  the  only  parameter  requiring 
time  discounting)  we  may  derive  X  ,  the  average  net  loss  per  unit  time,  as 
in  section  ?  and  write: 

6(k  4  P1  +  t.  +  b  ) 

X(R,  A)  -  d  -  M  ♦  d.  (R  -6  +  A  +  1)  ♦ _ _ _ _  4  £  (p.,  ♦  t0) 

*  ?  A  ' 

£  bn  r°° 

_  F(R)  >tcj  (x  -  R)  dF(x), 

A  A  R 


4  £  p-3a  - 


having  substituted  A  for  0.  Thus,  we  obtain 

l 


£(B'+b) 

(24)  X(R»A )  *  A*  +  - - -  + 


1a  +  d  R  -  _2  F(R) 
2  4  A 


„OQ 

+  6c  )  (x  -  R)  dF  (x), 
A 


where  A'  =  dj  -  M  +  f(p2  +  t?)  +  d4  f  > 


Bf  53  k  +  +  ^1 

dt2  “  d4  +  ?fp3 * 

Equation  (23)  is  of  the  same  form  as  equation  (5),  and  the  solution  for 
optimal  R  and  A  follows  as  before. 


p- 


PART  II  -  THE  SYSTEM 

9.  SYSTEM  PROCUREMENT  POLICIES 

In  this  part  of  the  paper,  several  concepts  of  optimization  on  a  system¬ 
wide  le^el  will  be  introduced,  where  a  system  is  defined. as  consisting  of 
several  bases  and  depots.  These  formulations  apply  to  both  non-recoverable 
and  recoverable  items,  and  will  be  set  forth  separately  for  life-of-type 
procurement,  periodic  procurement,  and  open-contract  procurement. 

Before  we  begin,  let  us  further  describe  these  three  kinds  of  pro¬ 
curement  policies.  Life-of-type  procurement  consists  of  the  purchase  at  one 
time  of  enough  of  a  spare  part  for  an  end  item,  such  as  an  airplane  or  a 
vehicle,  to  keep  the  end  item  in  repair  for  its  expected  remaining  life. 
Needless  to  say,  the  prediction  of  the  expected  life  of  an  end  item,  let 
alone  the  prediction  of  what  and  how  many  spare  parts  will  be  needed,  is 
an  almost  insurmountable  endeavor;  yet  such  a  procurement  policy  is  often 
warranted  in  view  of  the  enormous  retooling  and  setup  costs  to  make  parts 
for  an  end  item  no  longer  in  production.  This  is  particularly  true  with 
aircraft  which  have  a  notoriously  short  production  period.  Life-of-type 
procurement  may  also  be  used  for  very  low  cost  items  where  the  procurement 
costs  (cost  of  contracting,  etc.)  are  very  much  higher  than  the  cost  of 
the  items. 

Periodic  procurement  is  the  regular  purchase  of  supplies  to  satisfy 
demands  for  them  during  a  particular  fixed  period.  Thus,  w*=»  may  procure 
materiel  at  the  beginning  of  a  year  to  last  specifically  through  the  year 
and  supposedly  no  longer  then  the  year.  A  variant  of  this  policy  occurs 
’when  several  deliveries  are  specified  curing  the  period,  perhaps  at  equal 
time  intervals’ and  for  specified  amounts  at  each  delivery.  Periodic  pro- 
curement  and  its  variants  are  very  common  in  military  organizations  due 
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primarily  to  the  budgeting  and  authorization  of  funds  on  a  fiscal  year  basis. 

Open-contract  procurement,  in  its  simplest  font,  consists  of  procuring 
"on  demand",  or  when  the  stock  balance  reaches  a  certain  warning  level.  Thus, 
it  is  very  similar  to  the  reorder  point  concept  of  our  base  policy.  There 
are  many  variations  in  this  kind  of  procurement,  such  as  guaranteeing  minimum 
amounts  of  orders,  specifying  different  per  unit  prices  for  different  size 
orders,  etc.  This  type  of  procurement  is  very  convenient  for  avoiding  the 
high  costs  of  contracting,  but  is  more  or  less  limited  to  items  of  common 
supply,  or  items  where  retooling  and  setup  costs  are  relatively  low.  Also, 
in  military  organizations,  there  are  budgetary  and  authorization  restrictions 
on  the  use  of  this  procurement  policy. 

A  fourth  type  of  procurement  policy  is  possible;  one  which  is  not  now 
in  general  use.  This  is  the  case  where  th'e  procurement  periods  are  variable 
rather  than  being  fixed  as  in  the  case  of  periodic  procurement.  This  policy 
uses  a  reorder  level  to  trigger  procurement,  and  is  similar  in  concept  to  the 
base  formulation.  This  type  of  procurement  policy  is  discussed  with  the 
ca^e  of  open-contract  procurement,  since  the  formulation  is  identical.  It 
will  be  referred  to  as  "variable  period  procurement". 

The  decision  between  the  different  kinds  of  procurement  is  a  very  diffi¬ 
cult  problem  and  is  heavily  dependent  upon  the  particular  production  problems 
of  the  item.  We  do  not  attempt  to  solve  this  decision  problem;  we  assume  the 
type  of  procurement  to  have  been  somehow  already  established  for  each  item. 

We  are  probably  optimistic  in  saying  that  we  defer  thi3  problem  to  future  study. 

10.  THE  SHADOW-PRICE  CONCEPT 

In  the  equation  for  average  expected  cost  of  supply  at  the  base  per  unit 
of  time  (equation  (5)  in  section  2),  the  coefficienct  was  introduced  as 
that  part-  of  the  warehousing,  wear-and-tear  depreciation,  and  obsolescence 
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costs  which  vary  with  y,  the  amount  held  at  the  base*  This  definition  of 
d2  needs  to  be  reconsidered  when  the  base  is  treated  as  part  of  a  system 
of  bases  and  depots. 

First,  let  us  assume  that  the  item  incurs  the  same  warehousing,  wear- 
and-tear  depreciation,  and  obsolescence  costs  if  it  is  stocked  in  a  depot 
as  if  it  were  carried  at  the  base.  This  seems  reasonable,  particularly  in 
view  of  the  fact  that  the  obsolescence  cost  is  incurred  at  the  time  of  con¬ 
tracting  for  the  item  -  the  military  system  must  face  this  obsolescence  cost 
once  it  purchases  the  item,  regardless  of  where  the  item  is  stocked.  .  Let  us 
agree,  then,  not  to  charge  these  costs  to  any  of  the  bases  but  rather  to  a 
"system"  cost,  in  which  we  will  later  include  procurement  costs,  contracting 
costs,  and  other  costs  which  pertain  more  to  the  system  as  a  whole  rather 
than  to  any  component  of  the  system. 

On  the  other  hand,  if  it  does  not  cost  the  base  anything  to  hold  the 
items,  then  it  is  intuitively  evident  that  an  infinite  base  stockage  offers 

least  base  cost,  because  the  depletion  cost  becomes  zero,  and  the  paperwork, 

/ 

packaging,  and  transportation  costs  are  incurred  only  once.  This  is  also 

seen  by  inspection  of  equation  (8),  where  if  dp - ^  0,  A  - 7  OO. 

Therefore,  we  must  seek  other  factors  which  will  restrict  base  stockages. 

The  most  obvious  restricting  factor  on  base  stockages  is  the  fact  that 
there  is  only  a  finite  amount  of  the  stock  in  the  system  as  a  whole.  Let  us 
suppose  that  the  system  is  closed  in  the  sense  that  stocks  do  not  move  in  or 
out  of  the  system  after  procurement  has  been  made.  Then  it  seems  natural  that 
the  various  bases  should  "compete"  for  the  fixed  amount  of  stocks  available. 

So  that  they  may  compete  equitably,  we  impose  an  internal  "price"  for  the  item 
and  charge  this  price,  for  each  unit  of  time,  to  the  average  amount  stocked 
at  the  base  whenever  we  use  the  results  of  Part  I  in  calculating  base  levels. 
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After  base  levels  are  calculated,  we  do  not  include  this  factor  in  assessing 
base  costs,  since  it  is  in  fact  not  actually  incurred.  We  call  this  internal 
price  the  "shadow"  price  because  it  is  used  only  in  effecting  stockage  dis¬ 
tribution  ana  has  no  particular  relation  to  the  actual  procurement  price  for 
the  item. 

We  now  have  generally  identified  two  kinds  of  costs:  costs  associated 
with  the  bases  and  costs  associated  with  the  system  as  a  whole.  It  is  intui¬ 
tively  evident  that  the  sum  of  these  two  kinds  of  costs  is  minimized  by 
placing  all  the  available  amounts  at  the  bases.  Any  extra  amounts  at  the 
bases  certainly  reduce  the  base  costs,  while  the  costs  attributed  to  the 
system  remain  the  same  (after  procurement  has  been  made),  whether  the  items 
are  kept  at  the  depot  or  at  the  base.  Therefore,  we  are  in  a  position  of 
incurring  least  over-all  costs  by  having  no  depot  stocks.  .This  leads  us  to 
seek  costs  incurred  by  the  system  by  not  having  depot  stocks. 

If  we  put  all  the  stocks  at  the  bases,  then  a  particular  base,  depleting 
its  amount  on  hand  to  its  reorder  level,  can  obtain  a  resupply  only  from 
other  bases.  If  we  suppose  the  depots  to  be  the  original  source  of  the 
materiel  when  the  initial  distribution  to  the  bases  was  ma.de,  then  more  costs 
are  incurred  by  obtaining  resupply  from  other  bases  than  if  the  items  were 
left  at  the  depot  to  be  used  for  resupply.  In  other  words,  for  resupply  to 
a  base,  it  is  cheaper  to  leave  stocks  at  the  depot  for  the  resupply  than  to 
ship  them  to  a  base  and  then  later  ship  them  to  another  base  requiring  re¬ 
supply.  Any  shipment  from  one  base  to  another  is  an  excess  cost  in  this  re¬ 
spect.  We  call  the  costs  incurred  by  this  trans-shipment  of  goods  between 
bases,  the  "excess  transportation  costs",  even  though  they  include  extra 
packagin'',  an  l  other  associated  costs. 

Of  course  it  may  not  be  desirable  to  avoid  all  of  these  excess  trans¬ 
portation  costs.  If  we  hold  more  stocks  at  the  depots  we  cause  more  base 
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costs,  mainly  because  we  incur  higher  depletion  costs  at  the  bases.  On  the 
other  hand,  if  we  place  few  stocks  in  the  depots,  the  increased  base  stock- 
ages  provide  lower  base  costs  but  higher  excess  transportation  costs  caused 
by  trans-shipment  between  bases  after  depot  stocks  are  depleted.  There  is 
some  depot-to-base  stockage  ratio  which  should  minimize  expected  cost  to  the 
system  as  a  whole. 

Let  us  next  tie  these  concepts  together  and  express  them  in  symbolic 
forms.  We  begin  by  establishing  some  definitions  and  notations: 

(l)  Procurement  Period  «  length  of  time  between  two  successive  arrivals 

of  materiel  from  the  factory  after  procurement. 
This  definition  is  used  particularly  in  context 
with  the  periodic  procurement  policy. 

(?)  Lead  Time  **  length  of  time  between  the  calculation  of  system  require¬ 
ments  and  the  arrival  of  the  materiel  from  the  factory. 

(3)  Consumption  Period  55  lead  time  plus  the  subsequent  procurement  period. 

This  definition  is  used  particularly  in  context 
with  the  periodic  procurement  policy. 

(4)  t  «  number  of  units  of  time  in  the  procurement  period. 

(5)  d  e  system  shadow  price,  expressed  in  terms  of  the  same  unit  of  time 

as  the  procurement  period, 

(6)  d^  ®  system  shadow  price  ,  expressed  in  terms  of  the  routine  depot-to- 

base  pipeline  time  for  the  i-th  base. 

«  p^d,  where  is  the  routine  pipeline  time  for  the  i-th  base 

expressed  in  the  procurement  period  unit  of  time. 

(7)  Y  c  expected  total  base  stockages  at  the  end  of  the  consumption  period. 

(8)  Se  =>  expected  total  system  stockages  at  the  end  of  the  consumption 

period. 

(q)  s  *•  total  system  stocks  on  hand  at  the  time  of  reouirements  calcula- 
o 

tion,  or  at  the  beginning  of  the  consumption  period. 
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P  =  amount  procured. 

S  =  P  +  Sq  =  amount  which,  if  available  at  the  time  of  requirements 

calculation,  would  satisfy  system  demands  throughout  the 
consumption  period., 

J  =  excess  transportation  costs  incurred  during  the  procurement  period^ 
where  excess  transportation  costs  are  defined  as  transportation, 
packaging,  and  associated  costs  of  all  materiel  which  is  shipped 
from  one  base  to  another  base.  Transportation  costs  from  depot 
to  base  are  not  included  in  J, 

Now  let  us  look  at  the  base  costs.  First,  we  agreed  to  set  d2  to  zero 
in  the  base  cost  equation,  equation  (5),  in  order  not  to  charge  the  deprecia¬ 
tion  and  other  holding  costs.  We  also  do  not  include  any  factor  involving 
the  shadow  price,  because  this  is  used  only  in  calculating  A  and  R. 

The  actual  direct  costs  incurred  by  the  i-th  base  are  then  given  by 


(10) 

(11) 


(12) 


(25)  ^  (A^Rp  -  A!  + 


PBi 

~K~ 


£ib0i 


Cl  -  Fi^i)] 


+ 


where  A*  =  (p^  +  t^)  +  d-^,  and  the  other  terms  are  as  previous  defined. 

Here,  represents  only  the  fixed  set-up  costs  encountered  in  including 

the  item  in  base  inventory,  whereas  in  ecuation  (f>),  dj  also  included  the 
fixed  part  of  the  holding  costs,  assumed  to  be  a  linear  function.  In  eouation 
(?5),  \  and  Ri  represent  the  optimal  values  calculated  as  in  Part  I  from 


the  following  equations: 
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i  l  '-A"*1  '  c‘  -l  /"<*-%> 
(26)  Ri 


dFi(x) 


2  f  I 

A.  •  _  Bi+b01'  &-Fi(R10  +  J  (x  -  Ri)dPi (x)  l  . 

di  (  '  Ri 


Therefore,  both  A^  and  become  functions  of  d;  hence  If  in  equation 

(2 5)  is  a  function  of  d„  As  we  have  previously  noted,  X*(d)  increases  as 
A^  and  R,^  decrease  due  to  higher  expected  depletion  costs,.  This  implies 
that  ty(d)  increases  monotonically  as  d  increases*  For  the  i-th  base, 


X*  (d)  might  appear  as  follows: 


Figure  (8) 


The  function  may  or  may  not  be  concave  upward  as  shown. 

In  order  to  aggregate  the  values  for  X’  for  n  bases,'  we  must  first 

i 

convert,  them  to  a  comparable  unit  of  time.  It  must  be  remembered  that  Xf 

i 

(i  “  1  to  n)  is  the  expected  base  stockage  cost  per  unit  of  tine,  wnere  the 
tine  unit  is  the  routine  depot-to-base  pipeline  time  for  that  base.  There¬ 
fore,  we  divide  each  >f  by  p.. ,  the  depot-base  routine  pipeline  time 
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measured  in  terms  of  the  same  unit  of  time  as  the  procurement  period;  multi¬ 
ply  by  t,  the  number  of  time  units  in  the  procurement  period,  and  get  our 
total  expected  base  costs  for  the  procurement  period  to  be 

n  X! 

(27)  X  “t  Z  _  for  n  bases. 

3  i-i  Pi 

Of  course,  X  is  also  a  function  of  d  since  each  X*  is,  and  X 
s  is 

will  also  increase  monotonically  as  d  increases,  since  this  is  true  for 
each  X!  • 

The  excess  transportation  costs  <J  certainly  decrease  if  the  amount 
of  base  stocks  decrease,,  This  is  due  not  only  to  the  fact  that  fewer  base 
stocks  permit  less  trans-shipments  between  bases,  but  also  to  the  fact  that 
fewer  base  stocks  provide  more  depot  stocks  which  can  be  used  to  satisfy  the 
base  demands.  Since  we  have  previously  observed  that  base  stocks  decrease 
monotonically  as  d  increases,  we  conclude  that  J  decreases  monotonically 
as  d  increases.  The  excess  transportation  costs,  expressed  as  a  function 
of  d,  might  then  appear  as  follows: 

Figure  (9) 

J  I 


d 
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The  function  J  may  or  may  not  be  concave  upward  as  shown. 

The  excess  transportation  costs  also  depend  upon  the  amount  procured,  P. 
If  we  procure  a  large  amount,  then  we  increase  depot  stocks  and  decrease  the 
costs  J.  If  we  next  assume  all  other  costs  in  our  system,  namely,  the  depot 
and  procurement  costs,  to  be  functions  of  P  only  and  given  by  the  functions 
D (P) ,  then  we  can  write: 

(28)  L»(P,d)  *  D(P)  +  j(P,d)  +  X(d) 

s 

for  the  total  cost,  including  procurement,  of  operating  the  system  during  the 

procurement  period.  Theoretically,  values  for  P  and  d  may  be  found  which 

1 

minimize  this  total  system  cost  . 

The  factor  D(P)  in  equation  (28)  is  independent  of  base  stockage  dis¬ 
tributions  and  depends  entirely  on  the  procurement  costs  and  system  depletion 
costs.  Various  functions  for  D(P)  that  correspond  to  the  different  pro¬ 
curement  policies  will  be  discussed  in  sections  13  and  15. 

11.  THE  GENERAL  SYSTEM  SOLUTION  FOR  NON- RECOVERABLE  ITEMS 
In  seeking  values  for  P  and  d  which  minimize  the  total  system  costs 
expressed  by  equation  (28),  we  find  that  the  function  J(P,d)  for  excess 
transportation  costs  is  impossible  to  obtain.  However,  we  can  avoid  its  use 
by  examining  the  circumstances  which  create  these  costs;  namely,  the  ratio 
of  base  stocks  to  system  stocks.  Let  u3  consider  the  ratio 

(29)  E  -  Y(d), 

S 

e 

where  Y(d)  is  the  expected  base  stockage  at  the  end  of  the  procurement 
period  and  S  is  the  expected  system  stockage  at  the  end  of  the  procurement 
period.  In  eouation  (2d)  we  consider  E  to  be  a  fixed  number  less  than  or 

1 

Actually,  there  is  an  implicit  restraint  on  optimal  values  for  P  and  d; 
namely,  th  t  the  total  average  expected  base  stocks,  as  determined  b^  d,  do 
not  exceed  the  total  available  system  stocks,  as  determined  by  P.  .his 
constraint  is  expressed  by  setting  E  *  1  ir,  equation  (29,  . 
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equal  to  1.  This  equation  should  be  considered  as  a  means  of  " forcing”  the 
amount  procured  to  obtain  a  desired  base-system  stockage  configuration  at 
the  end  of  the  procurement  period.  The  constant  E  may  be  set  much  less 
than  1  to  provide  depot  stocks  for  war  reserves  or  "insurance"  or  to  conform 
with  other  command  policies.  Of  course  E  may  also  be  set  equal  to  1  to 
place  all  expected  system  stocks  at  the  bases  at  the  end  of  the  procurement 
period. 

We  may  now  express  the  system  costs  as  being 
(30)  L(P,d)  -  D(P)  ♦  Xs(d), 

subject  to  the  restriction  of  equation  (29).  As  will  be  subsequently  demon¬ 
strated,  equation  (29)  and  (30)  can  be  used  to  find  both  the  amount  of  pro¬ 
curement,  P  and  the  subsequent  stockage  distributions  for  the  various 
bases  and  depots.  The  remainder  of  this  section  will  be  devoted  to  obtaining 
eouations  from  which  the  procurement  amount  P  may  be  found. 

In  equation  (29),  we  may  express  Sg  as  equal  to  S  -i'  ,  or  P  +  S0  -X 
where  £  is  the  expected  consumption  during  the  consumption  period.  The 
value  is  the  system  stocks  on  hand  at  the  time  of  reauirements  calcula¬ 

tion,  and  hence  known.  The  value  £  is  assumed  to  be  a  constant  and  is 
evaluated  in  section  13  for  different  procurement  policies.  Since  enuation 
(29)  is  now  expressed  in  terms  of  d  and  P  we  may  find  the  value  of  P 
which  minimizes  L(P,d)  in  eouation  (30),  subject  to  the  restriction  of 
equation  (29),  by  the  use  of  a  Lagrangian  multiplier.  If  ju.  is  the  Lagrangian 
multiplier,  we  form  the  function 

P(f,d)  -  L(P,d)  ♦mDO)  -  E(P*S0-n] 

and  solve  the  equations 

(31)  a/1  -  QIHP,  -  liE  -  0 

ap  ap 
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aX(d)  «Y(d) 

IT  -  JL_  +  u  __ 

dd  ad  ^  ad 


simultaneously  with  equation  (29)  for  fi  ,  P,  and  d.  Eliminating  ja  in 
equations  (31),  we  thus  obtain  the  two  equations 

a  yd) 

E  _ _  +  6D(P)  6Y(d)  -  0 

ad  ap  ad 


(32) 


Y(d)  -  E(P  +  SQ  -X)  a  0 
to  be  simultaneously  solved  for  P. 

The  best  estimate  for  Y(d)  to  use  in  equation  (31)  is  given  by  the 

sum  over  all  bases  of  the  average  expected  base  stockages  as  given  by  equation 

(3).  Substituting  this  sum  in  equations  (32),  and  using  eauation  (?7)  for  Xs 

we  obtain  the  equations 

n  (  [R  (d),A  (d)] 

E-itg  11 

ial(p^  3d 

(33)  •  - 


6D(P) 

8P 


*3Ri(d)  3A.(d) 

___  - 


n 

v 


i“l 


R±(d) 


A. (d)  ♦  1 


E(P  *  So-j) 


•0. 


In  the  above  equations,  we  must  remember  that  R^(d)  and  A^(d)  are 
given  by  equations  (?6),  and  X?  is  given  by  eauation  (25).  The  function 
D(P)  is  derived  in  sections  13  and  15  for  various  procurement  policies. 

Equations  (33)  were  obtained  in  the  manner  shown  in  order  to  retain  the 
shadow-price  d  in  a  more  explicit  form.  Actually,  the  shadow-price  was 
used  only  as  a  heuristic  device  and  can  be  readily  eliminated  in  a  mathe¬ 
matical  solution.  Such  a  method  might  be  as  follows :  Eliminating  d^  from 
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equations  (26)  we  obtain: 


(33a)  <e .  (a . ,%)  -  =% +  L  O’  -  2Fi(Ri>  -  Aifi(RiQ 


"*oo  H  06  . 

2  f  (x  -  R^)dFi(x)  +  A^  JL  J  (x  “  %)dFi^x' 

4  1  m±  KJRi 


There  are 


actually  n  such  equations,  one  for  each  base,  each  with  different^ 


values  for  optimal  A  and 

We  are  now  interested  in  solving  the  equation 

n 

UP^.RJ  a  D(P)  4  t  1  - — 

1  1  i-1  Pi 

subject,  to  the  restrictions  of  the  equations  ^(a^)  “  0  and  the  equation 


tfWp)  '  ^  (Ri"£i + 


A,  +  1 

1  V  E(P  *  S  -J)  “  0 

2 


We  form  the  function 


Aa^R-.p)  -  UP.A.'V  +  l»i +^(VRi,F) 

where  p  and  p  '  are  Iagraifcgian  multipliers  (i  -  1  to  n). 

The  partial  derivatives  of  P  '  set  to  zero  are 


*lT  “  dD(P)-  p.E  “  0 
dP  dP 

ax,jL<K1»li)  ’’'VVV  .  0 

SC  *  -L  —Llf-  *  ■*!  -T1 —  11  0 

pt  eRj.  dHi 


dx’^(Rv»“p  (Ri’ \) 


+  ^  “  G 


1  dA. 

l 


m  t  x  r  ♦  u.  _ +  r  =  1 

a*-  — -  - rr - -  Hi  aa  2 

d4i  Pi  >  i 

Eliminating  the  Lagrangian  multipliers,  we  get  the  equations 
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ax»  a 4.  a<|  ax*  §pj  a  4 

(33b)  A  Et  _J; _ *  -  _ *  _ _  +  QD(P)  1  1  -  1 


Pi  ®i 


ap  7  0r.  av 

i  i 


which  must  be  solved  simultaneously  with  the  equations  <£'^(A^,P^)  “  0  and 
ip(Ai,Ri#P)  3  0  to  obtain  values  for  P,A^  and  PL. 

r)X  f  f  /7  QO 

Since  3.  -  -  h  B,  +  b  ,  Q  -  F.  (R. )]  +  0,  /  (x  -  R,  )dF,(x) 

3&  &2  1  01  11  'yR.  1  1 

1  1  U  1  -* 

and  -  -  JiCboifi(Ri)  ’  Ci  JL  f  (x  -  %)  dF.(xfl, 

ad  r—  ^  oi  l  l  i  «-  vy  1  i 


we  see  that 


WV 


SR.  -Rl 


ox;  at  ax; 

1  +  1  1  «  0 


FT  aS7  e,  Br7 

i  i  i  i 


ox!  0  au 

i  -  2  _ l  =  0. 

OR.  0A.  , 

Substituting  aij  lp  1  i  in  equation  (33t>)?  we  get 

9A,  2  3R, 

i  i 

~1^:“!jilrEt^i+  *P(F)  "]  «  0. 

_  2  aR±  aAi  J  L  Pi  ®i  ap 

It  can  be  shown  that 

1  H  -  X  CCifiOlj )  -  VPVl 

2  OR.  SA.  4  oi  1  i 

i  i 

Since  this  term  cannot  be  zero  for  all  R,  we  get 


at  i  +  ap(p)  »  0. 

p.  m±  ap 

The  system  solution  is  then  expressed  by  the  following  set  of  equations: 
Et  flXi  +  OD(P)  -  0 


(33c) 


Pi  3Ri 


8Xi  -2'9Xi 


S3 


0. 
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Z  (Ft,  *jL)  -  EP  +  2  (1  -  £ . )  +  E(<-  S0)  -  0 

i-1  2  1-1  2 

If  we  now  wish  to  re-introduce  the  shadow  price  d  we  observe,  from 


equation  (6),  that 


fb  .  f .  (R. )  -  (X  a 

*-  01  l'  1 


(x  -  R.  )dF.  (xQ  » 


*i  Ri 


1'  1 


Substituting  this  into  our  result  for _ i  we  see  that 

;»rj 

ax;  i 


which  must  hold  for  all  bases.  Therefore,  the  first  of  equations  (33c)  becomes 


1  »  i  ap(p)  m  d. 
pi  Et  a  p 

1  #  ax; 

The  first  of  equations  (26),  when  combined  with  our  result  for  shows  that 


so  that  the  second  of  equations  (33c)  is  no  longer  independent.  However,  in 
the  third  of  eouations  (33c.),  A  and  R.  may  be  expressed  as  functions  of 
d^  by  equations  (?6).  Therefore,  we  obtain  the  eouations 


(3  id) 


X  2D[P] 

Et  a  P 


n  A..(u  )  n 

2  (R,(J,)  *  ‘  1  )  -  E?  *  Z  (1 

L»1  7  1-1  7 


-  tj)  -  E(£-  S0)  »  0 


which  may  oe  solved  simultaneously  for  d .  an  i  P,  where  R^(d.)  and 
A«  (d-  )  are  given  by  equations  (T  ). 
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Equations  (33d)  represent  an  easier  set  of  eouations  than  equations  (33) 
for  finding  the  optimal  procurement  amount,  P. .  A  computational  method  might 
be  to  assume  a  value  for  P,  find  the  corresponding  values  for  d^  from  the 
first  of  equations  (33d),  and  calculate  optimal  A^  and  for  all  bases, 

using  this  value  for  d^  in  eouations  (?o).  When  these  optimal  values  are 
placed  into  the  second  of  equations  (33d),  together  with  the  assumed  value 
for  P,  a  result  is  obtained  which  may  or  may  not  be  zero.  If  different 
from  zero,  the  procedure  is  applied  for  different  values  of  P  until  a  P# 
is  obtained  such  that  the  result  is  nearest  zero.  This  P*  represents  the 
optimal  procurement  amount. 

Eouations  (33d)  provide  a  value  for  the  procurement  amount  P,  which  is 
"optimal"  in  the  sense  that,  given  E,  it  provides  least  expected  system 
costs.  Of  course,  there  nay  be  some  particular  E  which  yields  lower  costs 
than  any  other  value  for  E.  On  the  other  hand,  we  recognize  the  usefulness 
of  assigning  E  to  conform  with  military  policies.  In  a  sense,  then,  we 
obtain  a  sub-optimization  for  the  system  reouirements  calculation. 

12.  BASE  STOCKAGE  DISTRIBUTIONS 

Having  determined  the  amount  to  be  procured,  we  now  look  at  the  problem 
of  distributing  the  materiel  among  the  bases  and  depots.  Of  the  several  ways 
in  which  this  problem  can  be  approached,  we  consider  a  method  of  recalcula¬ 
ting  base  stockare  levels  at  a  number  of  times  throughout  the  procurement 
period. 

Let  the  procurement  period  oe  iivided  into  m  intervals  by  the  points 
t-j,  tp,  ...  tm+-j ,  where  t^  is  the  beginning  of  the  procurement  period  and 
tm+|  the  end.  At  each  time  of  calculation,  t,  (i  =»  1  to  m),  we  establish 


the  ratio, 
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n  A.(d)  +  1 

2  [R.(d)  -  £  *  J _ 1 

(34)  E  -  ,1-1  _ _ i _ ? 

S.  -  J. 

1  1 

where  the  numerator  is  the  expected  base  stockage  at  the  end  of  the  procure¬ 
ment  period,  is  the  total  system  stocks  at  t^,  the  time  of  calculation, 

l ^  is  the  expected  consumption  during  the  rest  of  the  procurement  period, 
and  E  is  the  same  constant  used  in  the  previous  section.  Various  values 
for  the  constant  JV  will  be  given  in  section  13.  Eouation  (34)  niay  be 

solved  uniquely  for  d.  This  provides  values  for  R.  and  A  (j  3  1  to  n 

d  j 

for  n  bases)  and  determines  the  base  stockages.’ 

The  distribution  of  amounts  procured  to  the  various  depots  is  determined 
at  time  t^.  The  procurement  amount  is  distributed  so  that  the  resulting 
depot  stocks  are  in  the  same  ratio  as  the  expected  mean  demands  of  the  bases 
normally  served  by  the  respective  depots.  The  mean  demand,  of  course,  must 
be  determined  in  terms  of  some'  common  unit  of  time.  If  there  are  two  depots 
A  and  B,  for  example,  then  the  respective  stocks  and  Sg  after  distri¬ 
bution  of  the  procurement  amounts  should  satisfy  the  relation 


M 

A 


st 


S 

B 


where  and  Mg  are  the  expected  mean  demands  of  bases  served  by  depots 

A  and  B  respectively. 

Let  us  now  look  at  the  effect  of  the  stockage  distribution  method  des¬ 
cribed  above.  If  part  way  through  the  procurement  period,  the  demand  for 
the  item  has  been  greater  than  anticipated,  the  expected  system  stocks  at 
the  end  of  the  procurement  period,  S|  -  f  j,  will  decline.  Therefore,  in 
order  to  maintain  the  ratio  S,  the  shadow  price  d  will  increase.  Thus, 
a31  base  reorder  levels  and  stock  control  levels  will  decrease  and  a  uniform 
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stringency  will  be  imposed  on  all  bases.  Base  stocks  will  reach  the  reorder 
levels  later  and  when  they  do.  smaller  orders  will  be  sent.  By  thinning  out 
stocks  at  all  bases,  this  procedure  reduces  the  probability  of  having  to 
backorder  any  of  them. 

If,  on  the  other  hand,  the  demind  for  the  item  has  been  less  than  antici¬ 
pated,  the  expected  system  stocks  at  the  end  of  the  procurement  period  will 
rise,  the  shadow  price  will  fall,  and  all  reorder  levels  and  stock  control 
levels  will  rise.  Thus,  base  stocks  will  reach  the  reorder  levels  sooner 
and  larger  orders  will  be  sent.  In  this  way,  the  extra  amounts  will  be 
spread  throughout  the  system. 

13,  THE  PROCUREMENT  COST  FUNCTIONS  FOR  EXPENDABLE,  NOW-RECOVERABLE  ITEMS 

In  this  section  we  consider  different  procurement  cost  functions  which 
correspond  to  the  different  kinds  of  procurement  policies  for  non-re coverable 
items.  In  this  section,  we  will  express  the  cost  functions  in  terms  of  S 
rather  than  in  the  procurement  amount  P.  Since  S  and  P  differ  by  the 
constant,  S0,  this  becomes  just  a  matter  of  convenience. 

First,  we  will  define  some  of  the  terms  to  be  commonly  used  in  this 
section. 

(l)  K  =•  fixed  cost  of  procurement,  including  contracting  costs  and 
s 

all  procurement  costs  which  do  not  vary  with  S.  For  open 
contract  procurement,  this  cost  is  the  fixed  paperwork 
cost  associated  with  placing  an  order  on  the  factory. 

(?)  U(P)  *  unit  procurement  cost  of  the  items  on  a  delivered  basis. 

This  function  includes  factors'-  to  depot  packaging,  in¬ 
spection,  and  transportation  costs.  When  P  “  0,  U(0) 
is  the  initial  set-up  cost 3  of  manufacture. 

(3)  c  »  scrap  value  of  the  items  or  the  salvage  value,  whichever  is 


higher. 
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(4)  p  *  routine  pipeline  time,  for  open  contract  procurement.  It 

s 

is  defined  as  the  number  of  days  from  initiation  of  an  order 
on  the  factory  to  the  time  of  arrival  of  the  materiel  in  the 
depot.  It  differs  from  lead  time,  as  previously  defined,  in 
that  it  contains  no  contracting  time  and  usually  little  or 
no  set-up  time  at  the  factory. 

(5)  0  ■  expedited  pipeline  time  for  the  relevant  procurement.  Its 

definition  is  the  same  as  for  p_  but  uses  premium  comoiuni- 
cation  and  transportation. 

(6)  f  (x)  «  system  demand  probability  functions  where  different  values 

i 

of  i  (i  a  1,  2,  3)  refer  to  the  different  procurement 
policies.  Thus,  f  (x)  represents  the  demand  probability 
for  the  remaining  life  of  the  item  as  used  in  the  life-of- 
type  procurement  policy,  f^Cx)  represents  the  demand 
probability  over  the  consumption  period  as  used  in  periodic 
procurement,  and  f^(x)  represents  the  demand  probability 
per  Pe ,  the  routine  pipeline  time  from  the  manufacturer 
to  the  depot.  This  last  function  is  used  in  the  open  con¬ 
tract  procurement  policy. 

rx 

(7)  F^(x)  “ J  f^(t)dt  a  cumulative  system  demand  probability  function 

o 

(i  -  1,  2,  3.)  corresponding  to  the  three- 

procurement  policies. 

rop 

(8)  £j a  J  xf . (x)dx  m  expected  mean  demand  (i  •  1,  2,  3)  correspond- 

■'o 

ing  to  the  three  procurement  p‘olicie3. 

(9)  *dj  +  Qe.y  *  cost  of  holding  an  average  amount  y,  where  the  sub- 

script  i  (i*l,  2,  3)  refers  to  the  three  procurement 
policies.  Thus,  when  i  ■  1,  this  cost  represents  the 


P-6/,7 

-56- 


holding  cost  over  the  remaining  life  of  the  item. 

When  i  “  2,  it  is  the  holding  cost  over  the  con¬ 
sumption  period,  and  ’when  i  ■  3*  it  is  the  holding 
cost  per  p  .  These  costs  include  interest,  ware- 
housing,  and  all  depreciation  costs  except  obso-  . 
lescence. 

(10)  b  +  b,  q  +  b„p„q  *  system  depletion  penalty,  or  the  cost  of 

insufficient  procurement.  In  this  cost, 

b  is  the  fixed  cost, if  an y,  which  is 
s  o 

independent  of  the  amount  of  depletion. 

The  cost  b.  is  the  cost  associated  with 
S  X 

the  amount  of  depletion  q,  but  not  dura¬ 
tion  of  depletion.  The  component  gb0pgq 
is  the  cost  associated  with  the  amount  and 
duration  of  depletion.  Here,  gb^  could 
be  the  cost  per  consumption  period  of  having 
an  extra  end  item  available  for  use.  The 
term  c  is  the  duration  of  depletion,  in 
terms  of  a  fraction  of  the  consumption  period, 
while  q  is  the  amount  of  depletion.  The 
various  components  of  the  depletion  penalty 
will  be  discussed  in  further  detail  or» 


A .  Life-Of-Type  Procurement 

The  procurement  cost  function,  D(S)  for  a  life-of-type  procurement 
policy  and  for  a  non-recoverable  item  is  given  oy 
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p  S 

(35)  D(S)  -  K  +  U(0)  f  Ql(S  -  SQ)  -  a]  I  (S  -  x)  dF  (x) 


+  d  + 
s  1 


)♦& 


/n\ 

vw 


+  Q  '  ?i(S)D  • 


The  second  tenn  Qj(S  -  3  )  -  cj  j  (3  -  x)  dF-^x),  in  the  right  mem- 
ber  of  equation  (35 )»  represents  the  expected  obsolescence  cost  of  procuring 
P  items#  In  this  term,  we  assess  obsolescence  of  an  item  procured  in  accord¬ 
ance  with  the  probability  of  that  item  not  being  used.  If  x  items  are  de¬ 
manded  during  the  life  of  the  item,  and  there  are  S  items  in  the  system 
initially,  then  S  -  x  items  are  not  demanded,  (x  <  S).  We  multiply  this 
number  of  items  by  the  unit  cost  of  the  item,  Qj(S  -  S  )  -  c~]  and  by  the 
probability  of  demanding  x  items  dF^(x) .  We  then  sum,  or  integrate,  for 
values  of  x  ranging  from  0  to  S  ,  where  the  symbol  S  means  integration 
to  S  but  not  including  S.  It  might  be  noticed  that  we  are  pricing  the 

amounts  S  which  are  on  hand  at  the  time  of  calculation  at  the  oame  unit 
o 

cost  as  the  amounts  procured. 


The  term 


in  equation  (35)  represents  the  average  expected 


? 

amount  on  hand  during  the  life  of  the  item.  Its  derivation  is  similar  to 
the  derivation  of  equation  (3).  This  term  is  multiplied  by  e.  the  com- 

S  x 

ponent  of  the  lifetime  holding  costs  associated  with  the  amount  held. 

The  last  term  of  equation  (35)  represents  the  depletion  penalty,  where 
U(0)  +  K  is  the  cost  of  depletion  and  1  -  F,  (S)  is  the  probability  of  de- 

S  #  *x 

pletion.  The  depletion  penalty  expressed  in  this  way  assumes  that  the  deple¬ 
tion  id  predictable  sufficiently  in  advance  to  obtain  more  materiel  from  the 
factory  before  actual  depletion  occurs.  Therefore,  we  limit  depletion  costs 
to  the  fixed  costs  of  procurement  plus  initial  set-up  costs  at  the  factory. 
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If  such  additional  procurement  is  made,  then  equation  (35)  is  used  again  to 
find  the  amount  of  procurement.  Of  couse,  the  function  F-^(x)  vdll  be 
different  in  this  case,  since  the  expected  remaining  lifetime  of  the  item 
has  changed. 

The  value  of  £  for  use  in  equation  (33)  is  £^.  To  find  the  value  of 
to  use  in  equation  (34)  we  suppose  that  is  given  for  i  **  1  to  in, 

where  £-^  is  the  expected  average  system  demands  during  the  i-th  calcula¬ 
tion  period,  and  there  are  m  such  periods  in  the  expected  life  of  the  item. 

n  m 

If  this  is  the  case,  £-^  “  4  "  Z  £^  an(*  -£•  *  Z  exPec^®d 

i»l  k°i 

demand  during  any  one  calculation  period  is  assumed  the  same  as  for  all 

£  n (m-i ) 

others,  then  £.  =  _ _ .  The  value  for  t  for  use  in  equations  (33) 

m 

is  the  expected  life  of  the  item,  expressed  in  days. 

B.  Periodic  Procurement 

The  procurement  cost  function  D(S)  for  a  periodic  procurement  policy 
and  for  a  non-re coverable  item  is  given  by 

S" 

r 

(36)  D(S)  ”  K  +  U(0)  +  IU(S  -  Sj  -  c  I  /  (S  -  x)  dF_  (x) 
s  ~  JQ  1 


+  d  +  e* 
s  ?  s  2 


)  *  b  ri  -  F,(S)  I 

p  so  £  ~ 


QO  ^QO 

+  A  I  (x  -  s)  dF  (x)  -  b„  J  (x  -  S)  (x  -  1  -  S) 

J-  O  *  a  c  O  r\  7.  7 


X  +  1. 


dFp(x). 


This  equation  is  similar  to  equation  (35)  except  for  the  last  three  terms 

which  represent  the  depletion  penalty.  The  term  sbJjL  -  F  (S)J  represents 

the  costs  of  depletion  which  are  invariant  with  the  quantity  of  depletion  and 

the  length  of  lepletion,  where  is  the  cost  and  jjL  -  F,  (S)^  tne  prob- 

r<*> 

ability  of  incurring  the  cost.  The  term  b_  J  (x  -  Sv  dF0(x)  represents 


P-647 

-59- 


the  costs  of  depletion  which  vary  with  the  amount  of  depletion  but  not  with 
the  length  of  depletion*  In  this  term  b  is  the  cost  per  depletion, 

3  X 

x  -  S  is  the. amount  of  depletion  and  dF^(x)  is  the  probability  of  demanding 
x  items  and  incurring  a  depletion  of  x  -  S  items*  The  integration  sums, 
for  values  of  x  from  S  to-  o o,  the  products-  of  amounts  of  depletion  by 

probability  of  incurring  that  amount  of  depletion. 

oo 

The  term  b  f  (x  -  S)  (x  +  1  -  S)  dF  (x)  represents  the  costs 
s  2  Js  2  (x+'i)  1 

of  depletion  which  vary  with  the  quantity  of  depletion  and  the  length  of 

depletion.  In  this  term,  b  is  the  per  unit  cost  of  such  depletion, 

s  < 

(x  -  3)  is  the  amount  of  depletion,  x  +  1  -  S  is  the  expected  duration 

¥  (ix  +  1) 

of  depletion,  and  dF?(x)  is  the  probability  of  depletion  by  such  an  amount. 
These  expected  costs  for  all  possible  depletion  amounts  are  summed  by  inte¬ 
grating  over  values  of  x  from  S  to  oo. 

The  depletion  penalty  defined  above  assumes  that  the  item  cannot  be  re¬ 
contracted  for  before  the  established  contracting  date  if  the  system  pre¬ 
maturely  runs  out  of  the  item.  However,  if  we  do  assume  that  additional 
quantities  may  be  contracted  for  before  the  regular  contracting  date,  then 
the  depletion  penalty  becomes  [K  +  U (0)3  [l  -  Fo(S)[],  as  was  the  case 
in  life-of-type  procurement.  This  assumes,  of  course,  that  the  system  is 
able  to  predict  the  impending  depletion  and  obtain  more  materiel  before 
actual  depletion  occurs. 

The  value  for  i  for  use  in  equations  (33)  is  As  in  the  case  of 

life-of-type  procurement,  tho  value  for  to  use  in  equation  (34)  is  given 

n 

by  /  «  EL.  where  is  the  exoected  mean  demand  during  the  k-th 

i  k.17<  -k 

calculation  period,  there  being  n  such  periods  in  tho  procurement  period. 

Tba  value  of  t  for  use  in  equations  (33)  is  the  procurement  period,  ex¬ 
pressed  in  days. 
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C.  Open  Contract  and  Variable  Period  Procurement 

For  the  open  contract  procurement  policy,  we  have  the  procurement  cost 
function  expressed  in  terns  of  two  variables  rather  than  the  single  variable 
as  in  the  other  policies*  These  two  variables  are  the  system  reorder  level 
R  and  the  system  reorder  amount  (amount  procured),  a  •  The  procurement 
cost  equation  is  very  similar  to  equation  (5)  for  the  base  costs.  The  same 
terminology  and  definitions  will  therefore  be  used  but  prefexing  s  as  an 
index  to  the  various  costs  to  distinguish  clearly  the  two  cases.  The  procure¬ 
ment  cost  function  for  open  contract  procurement  of  a  non-re coverable  item 
then  becomes 


£o(B  +  b  )  e0 

Dtt,.*)  +  -Li— _  +!i  A, 

A.  ? 


£,  b  EnC 

*  e„R  -  38  °  FjR„)  *  3  8  j 

3  3  s  3  3y  “  < 

As  Rs 


(x  -  Rs)  dF3(x) 


where  A  «  d  +  £  (  p  +  t  )  +  e  (1  -  £  ) 
s  s3  3  s  ?  s  2  s  3  j,  3 


B  *  K  +  p  +  t. 
s  s  si  s  1 


C  *  b.  +  b-  n 
s  8  1  3 


To  solve  for  optimal  R  and  a  we  must  develop  equations  to  replace 

s  s 

equations  (33 )«  An  eouation  analogous  to  equation  (30)  may  be  expressed  as 
follows : 


(38)  L(Rs>vd  )  >  D(RglAB)  ♦  X3(d  ). 


^his  equation  Is  subject  to  the  restriction  of  equation  (89),  where 


Se  ’  Rs  '  i  • 


1'sinq  tiie  metho-d  of  section  11,  we  form  the  function 


/  WRi’V  '  “(WW 
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and  solve  the  following  set  of  equations  simultaneously  for  Rg, and 


dZ  a  ap  -  o 

6As  dAs 


i/2  m  ML  -  “  0 

dR„  eft 

»  S 


“i 


+ 


P 


*»  0 


i/7  B  i  LLl..  + 

aAi  Pi  a-i 


Pi  ifL  +  u  ■  0 
a^:j  2 


^i(Ai,Ri)  *  0 


Using  the  results  of  section  11,  we  obtain  the  equations 


Pi  3D(RS^S) 

di  “ _ 

Et  3Rg 

(39) 

n 

£ 

i-1 


R^) 


\(di) 


n 

-  ER  *  E  (1  -(.)  *  E(  t-  S  ) 
s  i-1  ?  0 


0 


3D  =  0 

an 

'3 

which  may  be  solved  simultaneously  for  cf  f  Rg,  and  L\  Of  course,  R^(d^) 

and  ^(U^)  are  given  by  equations  (26). 

The  value  for  /  in  equations  (39)  is  £  .  The  value  for  t  in  equations 

(39)  is  p  .  It  i3  noticed  that  the  process  of  section  12  is  not  particularly 
s 

applicable  in  the  case  of  open  contract  procurement.  The  value  for  d  which 
determines  base  stockage  is  obtained  from  equations  (39),  and  does  not  change 
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unless  from  a  recalculation  of  equations  (39). 

Equation  (3?)  does  not  contain  obsolescence  costs  as  did  the  equations 

for  the  other  procurement  policies.  This  is  under  the  assumption  that  the 

open  contract  policy  would  not  be  used  when  the  item  approaches  the  end  of 

its  life*  However,  if  we  are  willing  to  recalculate  optimal  R  and  ^ 

S  8 

each  time  the  reorder  level  is  reached,  we  may  introduce  obsolescence  by 
adding  the  following  term  to  equation  (37): 


*3&(as)  -  <0 


R  +  t 
r  s  s 


(Rs  +AS  -  x)  dP1(x) 


which  is  as  previously  defined  except  for  dividing  by  *  9  to  express 

h 

the  cost  per  routine  pipeline  time,  p  , 

14.  THE  GENERAL  SYS TgM  SOLUTION  FOR  RECOVERABLE  ITg§ 

The  general  system  solution  for  recoverable  items  is  similar  in  concept 
to  the  solution  for  non-recoverable  items  presented  in  section  11.  Modifica¬ 
tions  must  be  made,  however,  to  allow  for  items  being  in  repair  and  hence  un¬ 
available  for  supply.  Also,  it  is  necessary  to  allow  for  gains  to  the  system 
from  attrittea  end  items  and  for  losses  due  to  condemnations. 

Before  introducing  the  system  cost  equation  for  recoverable  items,  further 
notation  must  be  developed.  This  notation  will  also  be  used  in  section  15. 

(l)  t  -  length  of  the  consumption  period,  expressed  as  a  number  of  units 
of  time.  The  unit  of  time  must  be  significantly  large;  in  the 
order  of  a  month.  If  the  unit  of  time  is  a  month,  then  there 
are  t  months  in  the  consumption  period. 

(?)  t*  *  the  remaining  lifetime  of  the  item,  measured  in  the  unit  of  time. 

(3)  ■  base  3tockag»  cost  for  the  i-th  base  (i=l  to  n)  during  the  j-th 

Pi 

unit  of  time  (j  *  1  to  t),  where  is  the  routine  pipeline 

time  of  the  i-th  bas^  exoressed  in  terms  of  the  unit  of  lime. 
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(4)  ■  expected  excess  transportation  costs  during  the  j-th  unit 

of  time  (j=l  to  t). 

(5)  r.  “  repair  cycle  ending  with  the  j-th  unit  of  time  (j»l  to  t)» 

.  t) 

Here,  r.  is  expressed  in  terras  of  the  unit  of  time# 

0 

(6)  f  (x)  83  system  demand  probability  function  over  repair  cycle  r.» 

4j  J 

x 

(7)  F,  .  (x)  •*  /  f,  ,  (t)dt  •»  cumulative  system  demand  probability 

4j  JQ  4j 

function  over  repair  cycle  r^. 

(8)  £  .  s  mean  expected  system  demands  during  the  j-th  unit  of  time 

4j 

(not  the  j-th  repair  cycle). 

(9)  w  *  average  wearout  rate,  or  the  ratio  of  items  condemned  to 

total  exchanged  items.  It  is  used  below  as  the  probability 
of  an  item  received  in  exchange  being  condemned. 

(10)  w  3  number  of  items  condemned;  a  random  variable. 

(11)  <|h(w)  «*  the  condemnation  probability  function  for  time  periods 

from  the  time  of  computation  up  to  but  not  including  re-^ 

pair  cycle  r.. 
poo 

“  j  <i»(x,w)  dF.(x),  where  <j,(x,w)  is  the  probability  of 
-'w  ^ 

condemning  w  items  when  x  items  are  demanded  (w  <  x), 

and  F.(x)  is  the  cumulative  demand  probability  function 
J 

for  time  periods  from  the  time  of  computation  up  to  but 
not  including  repair  cycle  r^.  It  can  be  shown  that  if 


x  takes  only  integral  values,  then  <|Ax,w)  “  (1  -  w) 
-w  .  1/ 


~\X-W 


wl  (x  -  w. 


1/  See  J.  V.  Uspensky,  INTRODUCTION  TO  MATHEMATICAL  PROBABILITY,  p  46. 
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(12)  ffij(w)  ='J  (j).(t)dt  **  cumulative  condemnation  probability 

function  for  time  periods  from  the 
time  of  computation  up  to  but  not 
op  including  repair  cycle  r^. 

(13)  <|)(w)  *  J  ^(x,w)  rdF1(x)  s  the  lifetime  condemnation  prob- 

w  r 

w  ability  function. 

r 

(14)  (D(w)  aJ  <|>  (t)dt  0  cumulative  lifetime  condemnation  prob- 

0 

ability  function. 

(15)  »  cumulative  condemnation  probability  function  for  the 

consumption  period. 

(16)  Sj  35  expected  total  number  of  items  (assumed  reparable)  obtained 

from  end  items  which  are  attritted  in  the  j-th  unit  of  time. 

It  does  not  include  condemned  items  from  attritted  end  items. 

Each  aj  is  assumed  constant. 

j~r, 

(17)  S*  ”  S  +  £  a,  D  expected  amount  of  serviceables,  reparables 

J  •  k-l  K 

and  condemned  items  in  the  system  at  the 
beginning  of  the  j-th  repair  cycle. 

(IB)  v(S)  =  average  expected  value  of  an  item  in  a  reparable  condition. 

This  value  may  be  given  as  the  procurement  unit  cost  ^(S), 
minus  average  cost  of  repairing  the  items  minus  average 
transportation  costs  from  base  to  depot.  The  value 
v(S)  -  c  is  always  positive  since  if  the  scrap  value  is 
greater  than  v,  the  item  is  condemned  rather  than  re¬ 
paired. 

(19)  p’  *  procurement  and  contracting  lead  time  measured  in  units  of 

3 

time. 


An  eouation  analogous  to  equation  (?tt)  ray  now  be  expressed  for  the  case 
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of  the  recoverable  type  item  as  follows: 


L(P)  *  D(P)  +  Z 


rS }  n5 j'W 

J  i  MjP.x.w)  dF  (x)  d  IB  (w) 


japf+1  W»0  X"0 

s 


t  fi 

+  I  j  M.(P,x«S»-w,w)  [1  -  F  (sH]  dffl.(w) 
w»0  J  J  4j  J  J 


j»p?+l  •wr.-Q 

s 


+  £  m,(p,x-o,w-s;)  Q  -  a.(s;Q 

j-pt+i  J  t  ~  i 


where 


n  >'  (P,x,w) 

H,(F,x,w)  -St  _  +  J.(P,x,w) 

^  i»l  ^ 


Equation  (40)  differs  from  equation  (28)  by*  the  addition  of  the  random 
variables,  w  and  x,  which  represent  the  amount  condemned  and  the  amount  in 
repair,  respectively.  The  function  M.  is  only  implicitly  a  -function  of 

J 

P,x,  and  w;  in  particular,  X*_,  is  explicitly  a  function  of  and  ^ .* 

Of  course,  equation  (40)  is  restricted  by  the  condition  that  total  expected 
base  stocks  do  not  exceed  total  expected  system  stocks.  This  restriction 
is  expressed  by  equation  (41)  below,  where  E  is  set  equal  to  1.  Different 
functions  for  D(P)  in  equation  (40)  are  developed  in  section  15  to  corres¬ 
pond  to  the  different  procurement  policies . 

In  equation  (40),  the  functions  Jj  are  impossible  to  obtain,  so  we 
again  resort  to  a  sub-optimisation  by  use  of  a  fixed  ratio  of  base  stocks  to 
system  stock.  Here,  however,  we  impose  this  ratio  for  each  unit  of  time  in 


the  consumption  period.  Therefore,  we  obtain  the  equations 

n  A.  .+1 

5&t1  ♦  _ -  qj 

(41)  e  -  e  -  i”1  iv _ 7 _ 


S»  -  x  -  w 

\) 
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By  a  method  analogous  to  that  of  section  11,  we  obtain  the  optimal 
procurement  amount 5  P,  as  a  solution  of  the  equation: 


p  S  t  p  S  ?  —  w 

J  3  J  3  d  .(P,x,w)  dP  (x)  d  ®  (w) 

,  _ 1J  4  J  J 


j«p*+l  w^O  x^O 


t 

♦  E 
j3’?  1*1 


r  S* 

J  J  d  .(P,x-S*-w,w)  Q- 
i  J  J 


F,  .(S!-w)]  dffl-j  (w) 
4j  J  J 


t  p. 

+  E  d  .  (P,x«0,w*S!)  Ql -0  .(S»)l  ■  _I 
j-pt+1  iJ  3  3  3  tE 


dP 


where  d  ,  (P,x,w)  are  given  as  solutions  of  the  eouations 
id 

n  f  A.  ,(d. ,)  “I  n 

E  R.  . (d  . )  +  1J  J-J  ■  t  E  (1  -  £  )  -  E.S!  +  E,(x+w)  «  0. 
i*8!  1J  1  J  ?  2  1J  J  J  d 

Of  course,  in  the  above  equations,  R,^(d  )  hnd  are  solutions 

of  equations  (8)  and  (9)  expressed  in  terms  of  the  i-th  base  and  j-th  period. 

j-r. 

Also*  it  should  be  remembered  that  S*  **  P  +  S  ♦  E  J  a,  . 

j  o  ,  -  k 

k«l 

The  above  results  are  used  only  to  find  the  procurement  amount  P,  and 
are  not  used  for  subsequent  base  stockage  distributions.  The  discussion  of 
section  12  still  applies  for  the  base  stockage  distributions,  except  that 
equation  (41 )  is  used  instead  of  eauation  (34)* 


15.  THE  PROCUREMENT  COST  FUNCTIONS  FOR  EXPENDABLE.  RECOVERABLE  ITEMS 


and  for  a  recoverable  item  is  given  by 


(43) 


D(S)  -  Ks  *  0(0)  +  Cv(S-S0)  -  o] 


t* 


t' 

(1  -  w  Z 


k«t i-r.  »+l 
t v 


'4k 


/  (St,  -  w)  dO(w) 


+  d  +  e 
8  1  si 


t»  t* 

E  ak~*  A  £4k  +  1 


S  +  k«l  k-1 


? 


(k  *  u(0)) 

“ 

t 

t  .1 

1  -  (D 

(s  ♦  s  a.  -  (1  -  w)  £  £,, 

9 

t 

s» 

k»l 

k«t-r  +1  J 

t  . 

b  Z 

A 

1  -  F  .  (S'-w) 

d  0  ,(w) 

8  ° 

J 

4j  j 

J 

♦  b  E 
3  1  J-P 1*1 


S» 

n  J 


.00 


(x-Sj+w)  dF^.(x)  d  fflj(w) 


Sf~w 

j 


s;  qo 

J  j  f1'  (x~S*+w)  (x+l-S»+w) 


♦  b  Z 

S  ‘  J-p»+l  o 
s 


J 

S*~w 

0 


dFy(x)  d®^(w) 


where  S^f  »  S  +  Z  a!c 


t» 

-  Z 


k-1 


,  £4k 


In  equation  (43)  the  second  term  in  the  right  hand  side  of  the  equation: 
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represents  the  cost  of  obsolescence.  In  this  term.  [v(S“So)  -  cf]  represents 
the  cost  of  one  item  becoming  obsolete.  By  using  Qv(S-Sq)  -  c^j  rather  than 
[U(M  )  -  cl  we  assume  that  the  system  is  able  to  predict  the  remaining  de¬ 
mand,  at  some  point  in  time  close  to  the  end  of  the  life  of  the  item,  sufficient¬ 
ly  well  to  halt  repair  on  items  that  are  not  going  to  be  issued.  Thus,  the 
system  saves  the  cost  of  repairing  and  transporting  the  reparable  which  is 
never  going  to  be  issued.  We  feel  that  if  [y(S--So)  -  cl  is  not  an  exact 
statement  of  the  loss  due  to  obsolescence  on  a  recoverable  item,  it  is  at 

any  rate  preferable  to  Qj(S-S  )  -  cl. 

t* 

The  term  (l-w)  Z  represents  the  number  of  items  that  must 

k®t  *-r,i+l 

become  obsolete.  It  is  the  number  expected  to  be  issued  in  the  last  repair 
cycle,  minus  those  that  will  be  condemned  in  that  time  period.  This  amount 
must  become  obsolete ' even  if  it  is  necessary  to  procure  more  of  the  item  to 
compensate. 

The  term  S  f  represents 'the  amount  of  stock  that  must  be  condemned  in 
order  to  avoid  all  obsolescence  on  the  item,  except  that  amount  which  cannot' 


be  avoided,  as. described  in  the  preceding  paragraph 
S.  . 


The  term 


{ 


(S^?~w)  d  ffi(w)  repregents  the  expected  number  of 


items  becoming,  obsolete.  The  term  S  -w  represents  the  number  of  items 
becoming  obsolete,  given  w  condemnations;  d(D  (w)  is  the  probability  of 


w  condemnations  in  the  lifetime  of  the  item,  and  the  product  (S  -w)  dffi(w) 

t  ^ 

represents  the  product  of  the  number  of  items  becoming  obsolete,  and  the 


probability  of  incurring  that  number  of  condemnations.  The  integral  sums 


these  products  frora  w  «  0  ,  in  which  case  the  whole  of  S  becomes  obsolete, 

to  5.  ,  ,  in  which  case  none  of  S  becomes  obsolete. 

t»  >  t* 


The  fifth  tern  in  the  ri$it  hand  side  of  the  equation: 
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e, 
s  1 


t» 

E  a. 


S  +  k-1 


t* 

-  w  E 
kEl 


f  *  1 


represents  the  holding  costs,  excluding  obsolescence,  where  the  bracketed 
part  is  the  average  expected  amount  of  the  item  on  hand  over  the  lifetime 
of  the  item. 

The  sixth  term: 


t  t 

(1-"}  =  +1 
k*l  k»t-r . +1 

t 

represents  the  expected  cost  of  having  to  procure  a  second  time.  The  term 
(K  +  U(0))  is  the  cost  of  having  to  procure  a  second  time; 

5 


(K  +  U(0)) 


1-33 


(1-w) 


t  _  t 

represents  the  probability  of  condemning  more  than  S  +  E  a,.  -  (l-w)  E  c 

k-l  k-t-rt+l  4k 

which  is  the  maximum  amount  of  stock  that  can  be  condemned  without  an  additional 
procurement  of  stock,  in  assessing  this  cost,  it  is  again  assumed  that  additional 
procurement  can  be  obtained  before  actual  depletion  of  serviceables  and  re par¬ 
ables  occurs. 

The  seventh  term: 


t 

b  E 
3  ° 


s« 

r  J 


j 


Cl  -  F, ,(s»  -  w)]  d  ffl(w) 

M-J  J  J 


represents  the  expected  system  depletion  costs  which  do  not  vary  with  either 
the  size  of  depletion  or  the  length  of  depletion,  resulting  from  a  temporary 


pile-UD  of  reparable  items  undergoing  repair.  The  term  b  represents  the 

s  o 
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cost  of  one  such  depletion;  the  term  [[l  -  F^  (3 *-w)[]  represents  the 
probability  of  exchanging  more  than  St-w  during  the  repair  cycle  r,; 

J  J 

the  integral  sums  for  all  values  of  w  from  0  to  S!j;  and  the  summation 

sums  the  depletion  probability  for  each  repair  cycle  r.  in  period  t,  Issues 

J 

(or  exchanges)  are  relevant  in  this  term,  since  each  item  in  the  system  at 
the  beginning  of  the  r  .-th  repair  cycle,  whether  serviceable  or  reparable 

J 

at  that  time,  is  available  for  issue  once,  and  only  once,  during  the  repair 
cycle.  If  the  item  is  serviceable  at  the  beginning  of  the  period,  it  is 
clearly  available  for  issue,  but  just  as  clearly,  the  item  which  it  replaces 
on  the  end  item,  whether  condemned  or  reparable,  cannot  be  repaired  before 
the  end  of  the  repair  cycle,  and  so  cannot  be  issued  during  the  period.  If 
the  item  is  a  reparable  at  the  beginning  of  the  repair  cycle,  it  can  be  re¬ 
paired  and  issued  before  the  end  of  the  repair  period,  and  will  be  available 
for  issue  once;  the  item,  it  replaces  cannot  be  made  serviceable  in  time  for 
issue  in  that  repair  period.  The  amount  S*-w  represents  the  stock  of 

J 

serviceables  and  reparables  on  hand  at  the  beginning  of  the  j-th  repair 
period. 

t 

The  summation  Z  represents  the  summation  of  independent  prob- 

j=,pUl 

abilities  of  depletion,  one  for  each  repair  cycle  r,.  Since  r  is  prob- 

J  J 

ably  larger  than  the  unit  of  time,  each  unit  of  time  may  be  in  several  re¬ 
pair  cycles.  The  probabilities  cf  depletion  in  each  repair  cycle,  despite 
the  overlap,  may  be  assumed  independent  if  the  unit  of  time  is  large  enough 
so  that  the  probability  of  depletion  in  one  repair  cycle,  in  -the  absence  of 
one  cf  the  units  of  time  therein,  is  negligible.  It  is  therefore  recommended 
that  three  units  of  time  be  the  maximum  in  ary  one  repair  cycle.  For  this 
reason,  the  month  has  been  suggested  as  the  unit  of  time.  The  lower  limit 


of  the  summation  is  o*+l  since  the  probability  of  depletions  occurring 

‘  s 
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prior  to  p?+l  cannot  be  affected  by  the  current  procurement  decision 
s 


The  eighth  term  in  equation  (43  )s 


t 

b.  E 
s  1 


S* 

'  J 


Q Q 


Sf-w 

J 


(x  -  S*  +  w)  dF^j(x)  d  {Dj (w) 


represents  the  expected  system  depletion  costs  which  vary  with  the  size  of 
depletion  but  not  with  the  duration  of  depletion,  resulting  from  a  temporary 
pile-up  of  reparable  .items  undergoing  repair,  here  b  is  the  cost  of 

S  JL 

the  depletion- per  unit  of  depletion?  (x-Sj+w)  is  the  size  of  the  depletion, 
given  a  Quantity  of  x  issues  in  repair  period  r.;  dF,  . (x)  is  the  prob- 
ability  of  X  issues  in  repair' period  r^;  th'e  second  integral  sums  the 
products  of  each  size  of  depletion  and  probability  of  that  size  of  depletion 
from  x  58  S*  -w  to  infinity;  and  the  first  integral  sums  the  products  of  the 
cost,  given  w,  and  the  probability  of  w  for  values  from  0  to  S*.  The 

J 

summation  is  the  same  as  the  summation  in  the  seventh  term. 

The  ninth  and  last  term  in  equation  (43) 


b  E 

s  ^  j«pt+l 
3 


Sj 


o  S*-w 


oo  (x-St+w)  (x+l-S!+w) 

_ a _ 1 _  dF,  (x)  d  CD. (x), 

?(x+i)  J 


represents  the  expected  system  depletion  costs  which  vary  with  both  the  size 
of  depletion  and  the  duration  of  depletion,  resulting  from  a  temporary  pile-up 
of  reparable  items  undergoing  repair.  The  term  b^  is  the  cost  per  unit  of 

S  •! 

depletion  of  a  depletion  lasting  throughout  the  repair  cycle;  the  term 


x+l-SS+w 

. _  is  the  expected  average  duration,  in  terms  of  a  fraction  of  the 

1  ?(x*i) 

re;>air  cycle,  of  a  depletion  of  an  amount  (x-S*+w);  the  term  (x-Sf.+  w)  is 

t>  J 

the  amount  of  depletion  given  x;  and  dF  (x)  is  the  probability  of  occurring 

4.1 

a  depletion  of  size  (x-S^w).  Tho  second  integral  sums  the  products  of  amount 

J 
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of  depletion,  duration  of  depletion,  and  probability  of  that  amount  of  de¬ 
pletion  for  all  x  from  St-w  to  infinity.  The  first  integral  and  the  sum- 

J 

mation  are  the  same  as  in  the  seventh  tern. 

B .  Periodic  Procurement 

The  procurement  cost  function  D(S),  for  a  periodic  procurement  policy, 
and  for  a  recoverable  item  is  given  by? 


Tne  difference  between  equation  (44)  and  equation  (43)  is  that  in  eoua- 
tion  (43),  t*,  the  remaining  lifetime  of  the  item,  and  t,  the  consumption 

period,  are  equal  while  they  are  not  ecual  in  equation  (44) •  In  this  equation 
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for  periodic  procurement,  we  assume  reprocurement  before  actual  depletion 
occurs.  However,  reprocurement  would  occur  only  because  of  condemnations 
and  not  because  of  pile-ups  in  repair.  Since  the  probability  of  condemna¬ 
tions  is  included  in  the  procurement  amount  calculation,  the  probability 
of  depletion  through  condemnations  becomes  relatively  small. 

C.  Open  Contract  and  Variable  Period  Procurement 

For  the  open  contract  procurement  policy,  and  variable  period  procure¬ 
ment  policy,  we  have  the  procurement  cost  function  expressed  in  terms  of  two 
variables  rather  than  the  single  variable  as  in  other  policies.  These  two 
variables  are  the  system  reorder  level  R  ,  and  the  system  reorder  amount 
Aa  . 

Before  introducing  the  cost  eouations  for  the  open  contract  and  variable 
period  procurement  policy  for  a  recoverable  item,  further  notation  must  be 
developed. 

(l)  p,  *  unit  of  time  *  period  from  the  time  system  stocks  (serviceables 

b 

plus  reparables)  reach  the  system  reorder  level  Rg  to  the  time 
the  materiel  is  received  from  the  factory. 

t* 

(?)  a  "  1  Z  a^  =  average  attrition  rate  per  unit  of  time  p  . 

t*  j-i  ■'  s 


(3)  £’  -  -  a 


net  expected  losses  per  unit  of  time  p  . 


average  system  procurement  period,  and  is  defined  as  the 
average  period  between  successive  deliveries  to  the  system. 


J-r,  _  j-r< 

(5)  S'!  -  R  *  t  +  E  J  a.  -  w  E  J  k , 
J  k=l  k=l 


(o)  Ail  other  terms  are  as  previously  defined. 

The  total  system  co3t  equation,  analogous  to  equation  (AO),  may  now  be 


expressed  as  follows: 


fs  +  1 
V 


(45)'  L(Rs>Aj)  -D(RS,4S)  ♦  £ 


^  VRs’Vx>  dF4,l(x) 

wrf 


&S  +  1 


MJ<V  V3^5  [1  -  FU(S«0 
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n 

where  M  ra  £  X:.(R  >A  »x)  +  d-s(RspAg,x^ 

J  i-1  1J  3  3  J 

Eouation  (45)  may  be  solved  for  optimal  Rg  and  Ag  by  setting  the 

partial  derivatives  with  respect  to  Rg  and  aq  equal  to  zero,  and  solving 

the  two  equations  simultaneously  for  Rs  and  Ag  ,  again  with  the  restraint 

that  expected  base  stocks  do  not  exceed  total  system  stocks. 

Again,  as  in  equation  (41 )  we  sub-optimize  by  replacing  the  Jj  functions 

with  the  restriction  that  R  and  a  satisfy  the  equations 

1J  i  j 


E  ■  E.  “  igl 

3 


l  [X, 


A  .  ,  +  l 

. .  ♦  ij  - 
ij 


We  may  now  write  the  procurement  cost  function  D(Rg/ig)  as  follows 


D(Rs,As)  =  _  +  gd3  *  se3  Rs  + 

S 


AS+1 

s  -  £* 


*  £'  (3P?  -  SV  *  ^  (sPl  *  Stl) 


S'  b  £* 

SO  V' 


A  i 
s  *’ 


zr«  r«r 

i  l  (x-5W 


s  ♦  1  oo 

£M^  &  r  (x-S^)  (x+l-S'j ) 


^  dF  (x). 
”  4j 


In  equation  (46)  the  first  four  terms  are  similar  to  the  same  terms  in 
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equation  (37)?  the  last  three  terms  are  similar  to  the  three  depletion  terms 
in  equation  (44),  -except  that  condemnations  are  treated  as  an  expected  mean 
value  rather  than  as  a  random  variable.  We  were  able  to  make  this  simiplifi- 
cation  in  equation  (46)  because  of  the  fact  that  the  system  reorder  level 
is  defined  in  terms  of  an  amount  of  uncondemned  items,  and  therefore,  con¬ 
strains  the  variation  in  the  system  stock  of  uncondemned  items  within  limits. 

Thus,  we  may  expect  S’J  to  attain  each  value  between  R  +  A  -  w£»  and 

J  s  g 

R  once  each  ordering  cycle,  and  may  expect,  by  averaging  over  a  large 
s 

number  of  cycles,  that  S"  will  exist  for  a  period  of  time  equal  to  that 

J 

which  would  be  obtained  by  using  an  expected  value  for  condemnations.  We 

(also  assume  that  is  shorter  than  the  repair  cycle,  so  that  any  demand 
which  occurs  after  Rg  is  reached  will  lead  to  a  reduction  in  the  stock  of 
serviceables  which  is  independent  of  whether  the  item  received  in  exchange 


is  condemned  or  reparable. 


Note  that  the  summation  in  the  depletion  terms  is  from  2  to  Qj*  1. 

* 

The  time  lag  of  one  period  represents  the  procurement  pipeline  p  ;  p  is 

s  s 

added  to  both  the  lower  limit  and  the  upper  limit  of  the  depletion  term 

saturations  since  the  R  and  to  be  calculated  cannot  affect  the  prob- 

s  s 

ability  of  depletion  during  the  first  period,  but  will  affect,  and  indeed 

A 

will  determine  the  probability  of  depletion  in  period  0  +  1  (or  s  +  1). 

S  £-}- 


If  we  are  willing  to  recompute  and  a  each  time  the  system  reorder 

^  s 

level,  R  is  reached,  we  may  add  an  obsolescence  term  to  equation  (46): 

s 

r  _  t»  o  i 

+  [y(S~S  )  -  cj  (1-w)  Z  +  f  (o-w)  dffi(w) 

L  k*t  ’-r.  t+l  -'o 

t» 

■  R  +  A  +  at*  -  t  £n( 

3  *  k«t*-rtVl 


where 


c 


P—647 
.  -76- 


In  computing  the  impact  of  the  procurement  requirements  calculations  on 
either  available  funds  or  budget  requirements,  the  calculation  is  carried 
forth  for  as  many  periods  t  or  9g  as  are  necessary  to  exhaust  either 
the  operating  or  the  budget  period.  In  the  life-of-type  calculations,  the 
procurement  amount  must  be  multiplied  by  the  unit  cost  to  obtain  the  funds 
or  budget  impact.  In  the  periodic  procurement  calculation,  the  procurement 
amount  must  be  calculated  for  as  many  procurement  periods  as  are  wholly  or 
partly  in  the  operating  or  budget  period.  In  each  calculation,  it  is 

necessary  to  estimate  the  balance  on  hand  at  the  time  of  the  next  calculation 

t-pf 

—  "  g  » 

This  amount  is  S’  ,  -  w  Z  £,  .  In  calculating  the  procurement  require- 
t_ps  k-1  4k  t-pt 

— ■ 1  x  cj  9 

ments  for  the  next  period,  SI  f  -  w  E  £  is  used  as  the  beginning 

ps  k»l 

balance  on  hand.  The  total  amount  to  be  procured  is  then  extended  by  the 
unit  cost  of  procurement  to  obtain  the  impact  on  funds  or  budget  requirements, 
For  the  open-contract  or  variable  period  calculations,  the  impact  on 
funds  or  budget  Requirements  is  Rg  +  n  Ag~  Sq,  where  n  is  the  number  of 
periods  9g  wholly  or  partly  in  the  operating  or  budget  peri  'd,  and  SQ  is 
the  current  balance  on  hand. 

If  there  is  a  fund  limit  within  which  all  procurement  requirements  must 
fit,  it  is  possible  to  sub-optimize  for  procurement  of  all  items  subject  to 
the  fund  limit.  Clearly,  the  fund  limit  must  be  considered  a  ceiling  only; 
that  is,  if  the  dollar  requirement  for  procurement  of  all  items  is  less  than 
the  fund  limit,  it  cannot  be  considered  as  either  optimal  or  sub-optimal  to 
raise  procurement  up  to  this  limit.  The  problem  before  us  ‘consists  of  the 
necessity  of  forcing  the  dollar  procurement  requirement  down  to  the  fund 
limit,  or,  in  other  words,  to  sub-optimize  procurement  under  constraint  of 


the  fund  limit. 
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For  this  purpose,  we  may  add  the  term  /1U(S)S  to  all  life-of-type  and 

periodic  procurement  equations;  and  the  term  fl  U(^  )  (R  +  A'  )  to  all 

s  s  8 

open  contract  and  variable  procurement  equations.  The  parameter  D  should 
be  held  at  zero  for  the  first  calculation  of  the  dollar  requirement  for  pro¬ 
curement,  If  the  total  dollar  procurement  requirement  is  greater  than  the 
fund  limit,  fl  may  then  be  raised  in  increments  until  the  dollar  procure¬ 
ment  requirement  is  equal  to  the  fund  limit. 

17.  REFINEMENTS  TO  THE  SYSTEM  SOLUTION 
A.  The  term  bQ  as  a  Function  of  Amount  of  Depletion 

In  all  previous  formulations  b0 ,  the  coefficient  of  the  cost  of 
system  depletion  which  varies  with  the  quantity  of  depletion  and  the  duration 

of  depletion,  has  been  assumed  to  be  constant.  We  may  wish  to  make  b  a 

s  2 

function  of  the  amount  of  depletion.  For  example,  the  military  organization 
may  be  able  to  prevent  any  loss  of  time  in  commission  of  the  applicable  end 
item,  despite  a  depletion,  through  the  device  of  "maintenance  cannibalization"; 
in  which  the  depleted  item  is  removed  from  an  end  item  entering  maintenance 
and  placed  on  an  end-item  which  only  needs  the  depleted  part  to  become  avail¬ 
able  for  use.  In  this  way,  an  end  item  out  of  commission  is  avoided  by  in¬ 
curring  the  cost  of  one  extra  removal  and  one  extra  installation  of  the  de¬ 
pleted  spare  item.  It  would  not  be  appropriate  to  use  this  cost  of  an  extra 
exchange  as  the  depletion  penalty  since  maintenance  cannibalization  is  prob¬ 
ably  limited;  it  may  suffice  to  prevent  end  items  out  of  commission  for  the 
first  ten  depletions,  but  beyond  that,  depletions  may  cause  end  items  to  be 

out  of  commission.  Therefore,  b  becomes  a  function  of  the  size  of  de- 

s  2 

pletion. 

We  may  introduce  this  modification  into  the  system  equations  merely  by 

moving  b.  inside  the  integral  as  a  function  of  the  amount  of  depletion 
s  2 
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■whenever  it  appears;  that  is,  in  equations  (36),  (37),  (43),  (44),  and  (46). 

For  example,  in  equation  (36),  the  last  term  would  become: 

» 

QO 

(  b  (x-S)  (x-S)(x+l-S)  dF  (x) 
s  s  1  2(x+l)  2 

Here,  bn(x-S)  is  the  cost,  which  varies  with  the  amount  and  duration 
s  l 

of  depletion,  of  (x-S)  depletions. 

The  same  modification  can  be  performed  on  h?  in  the  base  equations. 

B.  The  Obsolescence  Concept 

Our  treatment  of  obsolescence  in  the  system  equations  is  based  upon  the 
assumption  that  all  end  items  will  be  phased-out  as  whole  end  items;  that  is, 
not  lacking  the  spare  item.  Thus,  the  number  of  spare  items  condemned  repre¬ 
sents  a  fixed  minimum  requirement  for  the  spare  item.  The  logistics  system 
is  not  charged  the  procurement  cost  of  the  condemned  item  since  it  is  assumed 
not  to  be  within  its  discretion  whether  that  item  is  bought  or  not.  It  is, 
however,  charged  a  holding  cost  if  the  item  is  bought  too  early,  and  a  deple¬ 
tion  cost  if  the  item  is  bought  too  late.  The  assumption  will  be  completely 
valid  if,  for  example,  the  military  organization  wishes  to  mothball  the 
applicable  end  item  as  a  whole  end  item  at  the  time  of  phase-out;  or  to  pre- 
s  ent  the  end  item  to  a  friendly  nation,  again,  as  a  whole  end  item. 

If  the  military  organization  is  willing  to  phase  its  end  items  out  as 
incomplete  end  items,  we  may  substitute  the  term: 

U(S)S 

for  the  obsolescence  term  in  all  system  equations. 

C .  Variable  End  Item  Maintenance  Flow  Time 

The  terms  b  and  ,,’op,  whether  constant  coefficients,  or  variables , 
have  been  developed  on  the  assumption  that  the  maintenance  flow  time  of  the 
end  item  is  given.  If  b,.  and  gb^  are  interpreted  as  the  loss  of  utility 


F-4' 
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of  the  end  item,  they  must  be  modified  by  the  expectation  that  the  end  item 
requiring  the  part  may  have  been  out  of  commission  already  when  the  demand 
for  the  spare  item  arose.  The  end  item  may  have  been  out  of  commission  for 
another  part,  or,  more  likely,  may  have  been  out  of  commission  because  the 
need  arose  while  the  end  item  was  undergoing  overhaul.  If,  on  the  average, 
three  depletions  out  of  ten  lead  to  end  items  out  of  commission  and  seven 
out  of  ten  do  not, 'for  either  of  the  foregoing  reasons,  b  and  b  must 

A  3  <£ 

be  developed  as  .3  times  the  lost  utility  of  the  end  item,  rather  than  the 
whole  lost  utility  of  the  end  item.  If  the  maintenance  flow  time  of  the 
end  item  were  reduced,  as  a  matter  of  policy,  then  we  should  expect  that  a 
larger  per  cent  of  depletions  would  lead  to  end  items  out  of  commission. 

If  maintenance  cannibalization  is  used  to  measure  the  depletion  penalty, 
the  maintenance  .flow  time  of'  the  end  item  becomes  important  in  determining 
the  number  of  exchanges  in  the  spare  part  repair  period  of  time  in  which 
gbp  is  defined.  For  example,  if  the  spare  part  repair  period  of  time  is  30 
days,  and  the  end  item  maintenance  flow  time  is  10  days,  a  system  depletion 
of  the  spare  part  lasting  one  period  of  time  will  be  valued  at  the  cost  of 
three  exchanges  through  cannibalization.  If  the  maintenance  flow  time  f or 
the  end  item  is  then  cut  to  5  days,  the  same  system  depletion  of  the  spare 
part  will  be  valued  at  the  cost  of  six  exchanges. 

Thus,  b9  and  b  ,  and  also  L,  the  system  cost  of  logistics  support 
occasioned  by  one  spare  item,  become  functions  of  the  end  item  maintenance 
flow  time. 

Let  us  introduce  some  new  terminology  at  this  time : 

(1)  iRj  *  maintenance  flow  time  of  the  applicable  end  item  in  period  j. 

(2)  L^(m^)  «  lifetime  system  logistics  cost  of  spare  item  i  as  a  function 

of  m^. 
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(3)  t  ”  total  useful  life  of  the  end  item. 

(4)  M^(ra  )  55  lifetime  cost  of  a  pool  of  end  items  out  of  commission  for 

cl 

maintenance  as  a  function  of  m^. 

p 

(5)  Mam.)  =  cost  of  maintenance  on  the  end  item  in  period  j,  including 

direct  and  indirect  labor  and  capital,  as  a  function  of  m.» 

J 

(6)  «  total  lifetime  cost  of  maintenance,  logistics  support,  and  a  pool 

of  end  items  out  of  commission  for  both  maintenance  and  parts, 
over  the  lifetime  of  the  end  item* 

We  may  now  write  an  equation  for  as: 


n  1  t  p 

(47)  C,(m.)  -  t  L.  (ra.)  +  M  (m. )  +  l  M;(m  ) 

^  J  1  J  J  lei  J  J 


where  there  are  n  spare  items  applicable  to  the  end  item. 

We  may  then  take  t  partial  derivatives',  one  for  each  of  the  m.;  set 
them  equal  to  zero;  and  solve  simultaneously  for  the  m.» 

D.  The  Batching  Concept  in  Repair 

In  defining  the  repair  cycle  in  period  j  as  a  constant,  we  have  assumed 
that  all  reparable  items  are  scheduled  through  repair  as  fast  as  they  are  gen¬ 
erated  (although  the  repair  cycle  may  vary  between  periods).  The  military 
organization,  however,  may  be  faced  with  a  large  fixed  cost  per  batch  of 
items  scheduled  through  repair.  We  present  below  an  alternative  formulation 
of  the  system  cost  equation  which  may  be  used  for  determining  simultaneously 
the  requirements  from  procurement  and  the  optimal  batch  size  for  repair.  We 
shall  present'  the  case  of  ilfe-of-Hype  procurement  as  an  example. 

Let  us  .first  introduce  some  new  terminology: 

(l)  t*  »  the  lifetime  of  the  item  measured  in  major  periods  of  time 


(periods  in  the  order  of  a  year)  t  . 


P-6  L? 

-8X- 


(2) 

(3)  P 


the  major  period,  of  time  measured  in  minor  periods  of  time,  p  • 
the  minor  period  of  time  which  is  defined  as  the  period  from 
starting  the  batch  of  reparables  through  repair  to  their  emer¬ 
gence  as  serviceables. 


(4)  A^.  *  the  repair  quantity  or  batch  size  in  the  j-th  major  period;  (A^  -  1) 

,r)  5  «  s  +  ^  a +  fg  -  w  -  the  average  quantity  of  serviceables 

^ '  R  j  k3!  R  p  ' 

and  reparables  in  the  system  in  the 

j-th  major  period  (assumed  constant 

over  the  major  period),  who-e-  w,  as 

before,  represents  condemnations  and  is 

a  random  variable. 

/, s  .  „  „  «  the  repair  level,  or  quantity  of  serviceables  in  the 

;  Rj  ‘Rj 

system  at  the  time  that  repair  on  the  batch  is  oegun. 

(7)  f  (x)  -  the  demand  probability  function  over  the  minor  period  pR 

Rj 

in  the  j-th  major  period. 

(8)  FR.(x)  m  f  ^.Wdt  «  the  cumulative  demand  probability  function 

over  minor  period  pR  in  the  j-th  major  period. 

/Q\  2  «  ^  *  the  average  cycle  between  receipts  from  repair  in  the  j-th 

'  '  Rj  - 

Rj  major  period,  where  £  Rj  is  the  mean  of  fRj(x)* 

(10)  (D .  (w)  -  the  cumulative  condemnation  probability  function  for  periods 

prior  to  major  period  j,  and  for  half  of  major  period  j. 

(11)  p.  *  the  routine  depot -to-base  pipeline  for  base  i,  measured  m  units 

of  PR. 

n  v  .  ()  .  the  holding  cost  over  the  major  period,  where  y  is 

iR1  jrou' 


the  average  amount  of  the  item  in  the  system  in  the 
major  period 
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(13)  rp  +  ry(q)  =  the  cost ' of  repairing  a  batch  of  items,  where  a  is 

the  number  of  items  in  the  batch 

(14)  Jj (S^j,x,w)  18  exce<3s  transportation  costs  defined  over  the  minor 

period,  p^. 

An  equ&tion  analogous  to  equation  (40)  may  now  be  expressed  as  follows: 


t»  '  t  £ 

L(S,abJ  *  d(s,a  )  +  2  ■»  RJ 

Rj  RJ  J"P|+1  arj 


|M>  M1(SS-1,X,W)  Fg j  (x)dx 
\r*  0  Lu  ^ 


iW^", 


VW**w) 


[l  -  F(x'0  dx 


SRJ'VW 


(x-S  ,+^_  +w)  / 

■  r-1  Rj  .  MjCS^.jX-S^-Zi^-w^)  dF(x)  >■  d(D(w), 


where  H^(Sn,,x,w)  «  I  ^  ^  Rj*  *  ^  +  J4(S„.txtw), 

3  i-1  p.  J 

1 


Of  couse,  this  expression  is  subject  to  the  restriction 
n  A.  +1 

2  (Rjj  +  Jl_-  ) 

E  -  E,  -  i-1  XJ  P  J 

}  "  ""  — '  "  m - - - - - -  * 

SRj  -  x  '  * 

The  system  procurement  function,  D(S,A^)  may  he  expressed  as  follows: 


t*  S_. 

d(s,a  )  =  e  rR; 

ft.i  j=p*+i  j  v^c 

j 

+  ^  R.l^rc  rr  +  ^RjSn 

"rj  ij 


dRl  +  dR?  +  'Ttj  t,n  rv 


sbo&-FRj(3Rj-Apj-W)] 


♦  b 

s  1  WH 


(x-VVW)  dFaj(x) 
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'  -T. 


§ffl(w) 


i  I 

i 

I 


+  Tk  +  u(o)l 


03  (s 


t» 

+  £  aj, 

k=.l 


t* 

(1-w)  Z 


k”t 1 -r,  +1 
t- t 


,  4k  * 


There  are  two  rather  questionable  assumptions  implicit  in  these  equations. 
First,  it  is  assumed  that  the  wearout  rate  is  low  enough  for  S-^  to  be  a 
reasonably  close  approximation  to  the  stock  of  serviceable s  and  reparables 
throughout  major  period  j.  Second,  it  is  assumed  that  the  probability  of 
x  >  is  negligible.  The  formulation  could,  however,  be  refined,  and  the 
first  of  those  assumptions  relieved  somewhat,  by  reducing,  through  an  iterative 


process,  the  major  per* -4  tv  until  it  approaches  the  largest  0. 
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